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The ordinary surface magnetic phase transition is studied for the exactly solvable anisotropic spherical 
model (ASM) , which is the limit D — > oo of the D-component uniaxially anisotropic classical vector 
model. The bulk limit of the ASM is similar to that of the spherical model, apart from the role 
of the anisotropy stabilizing ordering for low lattice dimensionalities, d < 2, at finite temperatures. 
The correlation functions and the energy density profile in the semi-infinite ASM are calculated 
analytically and numerically for T > Tc and 1 < d < oo. Since the lattice dimensionalities d = 1, 2, 3, 
and 4 are special, a continuous spatial dimensionality d' = d — 1 has been introduced for dimensions 
parallel to the surface. However, preserving a discrete layer structure perpendicular to the surface 
avoids unphysical surface singularities and allows numerical solitions that reveal significant short- 
range features near the surface. The results obtained generalize the isotropic-criticality results for 
2 < d < 4 of Bray and Moore [Phys. Rev. Lett. 38, 735 (1977); J. Phys. A 10, 1927 (1977)]. 
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I. INTRODUCTION 

Magnetic ordering in semi-infinite and film geometries 
is an old problem currently receiving increasing atten- 
tion because of enormous progress in fabrication of mag- 
netic structures on the atomic scale. Theoretical methods 
using the mean field approximation (MFA) or the phe- 
nomenological Landau theory ||l|,||,^,|[^] provided classifi- 
cation and description of the qualitative features of differ- 
ent types of surface phase transitions. High-temperature 
series expansions [^| and Monte Carlo simulations 
as well as the scaling analysis and the e expansion 
|p|,PHlO|], shed light on the details of the surface critical 
behavior. A general review of these approaches can be 
found in Refs. |jll|,|l^ . Examples of recent work in film ge- 
ometry based on e expansion are Refs. [^|4j. A special 
case is the confined two-dimensional Ising {S ~k) model, 
for which exact solutions have been found [ |r5| , p^ , p7|JT^ . 

The ordinary surface phase transition of the semi- 
infinite ferromagnet occurs at the bulk critical tempera- 
ture Tc- It is characterized by a number of surface crit- 
ical exponents the definition of which can be found in 
Ref. Ill . In particular, the susceptibilities at the sur- 
face with respect to the fields applied either in the bulk 
or at the surface are described by the exponents 71 and 
711, respectively, which in the MFA are given by 71 = ^ 
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and 711 — — ^ (no divergence), in contrast to the bulk 
exponent 7=1. In thin films, which are more important 
for applications and more interesting for the experiment, 
there are additional effects, such as the lowering of Tc in 
comparison to its bulk value and the crossover between 
the three- and two-dimensional behavior as a function of 
the film thickness L |9|,|o|j2|]. 

The main body of the theoretical work on surface phase 
transitions is being done, with an exclusion of the 2d 
Ising model, starting from the field-theoretical contin- 
uous Hamiltonians or free energies. Such an approach 
has proven to be very useful for establishing the univer- 
sality classes and critical laws but, on the other hand, 
the (nonuniversal) absolute values of observables, such as 
critical amplitudes, remain undetermined. In addition to 
the restrictions of the field-theoretical methods that are 
well known from the bulk physics, there are more specific 
problems related to the role of the lattice discreteness in 
confined geometries. One can question how the contin- 
uous approach can be applied to thin films consisting of 
a mesoscopic number of layers. A similar question can 
be addressed to the semi-infinite ferromagnets as well — 
does the continuous approximation apply in the region 
near the boundary, n ^ 1, where n = 1,2,... is the 
layer number? The MFA, or the Landau theory, gives 
a positive answer to this question near criticality, where 
the correlation length is much larger than the lattice 
spacing ag and the order parameter — the magnetization 
— cannot change at distances smaller than ^c- This is, 
however, not the case if one goes beyond the MFA and 
considers spin-spin correlation functions. If the tempera- 
ture is high enough, or the system is classical, spin waves 
with the wave vectors up to the edge of the Brillouin 
zone k ~ fcmax = ""/ao are excited. This means that cor- 
relation functions comprise oq as the length parameter, 
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additionally to ^c, and thus there can be inhomogeneities 
near the surface of a ferromagnet on the scale of several 
atomic layers, even near criticality. This boundary re- 
gion, n ~ 1, is that which can be locally probed in ex- 
periments, and here the continuous approximation may 
become, at least quantitatively, wrong. 

Spin-wave effects in weakly anisotropic systems dras- 
tically change their behavior in comparison to the MFA 
predictions. Magnetic models with continuous symmetry 
in low dimensions, c? < 2, cannot order at finite temper- 
atures, and for models with d > 2 the correlation length 
is infinite in the whole region below Tc- These effects 
are not less important than the critical coupling of fluc- 
tuations giving rise to nonclassical critical indices. In 
Heisenberg systems the linear spin- wave theory satisfac- 
torily describes the above mentioned effects well below 
Tc but breaks down at elevated temperatures. There is, 
however, a model where a kind of spin-wave description is 
valid in the whole temperature range, whereas the critical 
fluctuation coupling vanishes. This is the D-component 
classical spin vector model proposed by Stanley [ ^J23| 
in the limit D ^ oo. Stanley has shown |24| that in 
this limit the partition function of a homogeneous ferro- 
magnet coincides with that of the spherical model (SM) . 
The latter was advanced by Berlin and Kac j2^ as an 
exactly solvable substitute for the Ising model. The for- 
malism contains, however, the spin- wave integral over the 
Brillouin zone and describes rather the properties of the 
isotropic Heisenberg model. The spherical model in its 
traditional form was extended to inhomogeneous systems 
by Barber and Fisher [ p6| , who found a nonmonotonic 
dependence Tc{L) for thin films. This unexpected fea- 
ture was attributed to the failure of the global spin con- 
straint in the inhomogeneous case. Later, an improved 
version of the SM was proposed, which used spin con- 
straints on each lattice site |2^, and it was shown that 
this version is equivalent to the D — oo model in the 
inhomogeneous case. Application of the SM with con- 
straints in each layer [ p8[ yielded a reasonable monotonic 
dependence of the numerically calculated Tc{L) for the 
films in four dimensions. Because of the complexity of 
models of this type, most researchers still use the more 
convenient global-constraint SM in confined geometries 
(see, e.g., Refs. p9| , ^^ ). Very recently, a compromise 
model was proposed||3^ which uses a constraint for the 
spins at the surface in addition to that for the bulk ones. 
The properties of such a model are closer, in a sense, to 
those of the SM with the local spin constraint and, hence, 
to those of the D = oo, or 0(oo), model. 

A remarkable property of the Stanley model is that it 
can be easily extended to the anisotropic case, whereas 
the traditional SM cannot. This is very important be- 
cause the anisotropy breaking a continuous symmetry 
plays a crucial role in low-dimensional systems, where 
it stabilizes ordering at finite temperatures. The limit 
_D — > oo of the uniaxial £)-component vector model deter- 
mines the so-called anisotropic spherical model (ASM), 
which is described in the inhomogeneous case by a closed 
set of equations for the variables on the lattice sites ob- 



tained in Ref. [g3|. The ASM was applied in Ref. |3|] 
to study the dependence Tc{L) for ferromagnetic films in 
three dimensions. Here, for all fixed values of the film 
thickness L one has Tc ^ in the isotropic limit because 
of the two-dimensionality of the film, and such a behav- 
ior is pertinent for any Heisenberg film with finite D as 
well. Thus, allowing for the anisotropy is relevant here, 
and it cannot be done within the traditional SM. An in- 
teresting feature of the solution obtained in Ref. [Q is 
the role of the correlation length of the transverse spin 
components, ^ca, in the crossover from the three- to two- 
dimensional behavior of the film, which takes place for 
L < ^ca- Note that in the finite-size scaling analysis (see, 
e.g., Ref. [pl| ) only the diverging longitudinal correlation 
length fez is used, whereas the noncritical transverse cor- 
relation length is disregarded as an irrelevant variable. 

Another application of the ASM is to the temperature- 
driven phase transition between the Bloch and linear 
(Ising-like) domain walls in uniaxial ferromagnets at 
some Tb < Tc This phase transition was studied 

within the framework of the MFA and the Landau the- 
ory in Refs. ]35|,|3^ ,^ , ^ , ^ and with a field-theoretical 
method in Ref. |p^(]|| 7^ A low minimum of the domain- 
wall mobility at Tb predicted in Refs. |4^,^^ was used 
to identify the domain-wall phase transition in the dy- 
namic susceptibility experiments on Ba and Sr hexafer- 
rites 1^,0 . It is clear that the anisotropy is an impor- 
tant characteristic of the model, giving rise to the very 
existence of domain walls of finite width separating the 
"up" and "down" domains. For this reason, the attempts 
to describe domain walls with the traditional SM in both 
versions with and without the global spin con- 
straint could not yield relevant results. 

In the recent work of Ref. it was shown that the 
ASM and SM are not equivalent, even in the isotropic 
homogeneous case, if the longitudinal correlation func- 
tion (CF) SzzO^) below Tc is involved. Whereas in the 
traditional SM the CF is proportional to 1/fc^ at small 
wave vectors, the behavior of 6*22 (k) in the ASM shows 
a more complicated behavior which is governed by the 
spin-wave effects and is sensitive to the dimensionality. 
In three dimensions (k) oc 1/fc, which is familiar from 
the linear spin-wave theory. The above law holds for 
k^Krn oc Tc — T , i.e., there is a critical length scale 
f„i = 1/Km in the theory. The length f^ is analogous 
to the "bare", or the mean-field, correlation length be- 
low Tc, whereas the true longitudinal correlation length 
fez remains infinite in the isotropic model below Tc. The 
former is responsible for the crossover of the real-space 
CF from Szz{r) (x 1/r for r<f„ to Szz{r) oc l/r^ for 



The ASM equations of Ref. |33|, as well as those for 
the SM without the global spin constraint |2|] , are rather 
complicated, strongly nonlinear equations for the vari- 
ables on a lattice. In the latter case some researchers 
termed them analytically intractable. Nevertheless, for 
the weakly anisotropic ASM in the domain- wall geometry 
it was possible to guess the solution |Q, which yielded 
an example of the phase transition of an interface that 
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is analytically tractable beyond the MFA . Surfaces with 
free boundary conditions make the problem much more 
complicated. In Ref. only the most important and 
partially rough asymptotes for the Tc shift in films could 
be obtained, and numerical calculations have not been 
performed. The aim of this work is to investigate the in- 
fluence of surfaces on magnetic ordering in more detail for 
the ordinary phase transition in the semi-infinite ASM in 
the temperature range T > Tc. As we shall see, analyt- 
ical solutions are available in the dimensionality ranges 
1 < d < 2 and d > 4, as well as for 2 < d < 4, both 
at and away from the isotropic criticality. In addition, 
the problem will be solved numerically in all the relevant 
cases. 

The same problem presented here was addressed in the 
seminal work by Bray and Moore , who considered a 
field theory with the volume and surface Hamiltonian 



densities of the type (V(/)) 



ucj)"^ and ccf)^ , re- 



spectively, for the n-component vector order parameter 
(j) in the limit n oo. Since this model is isotropic, 
the range 1 < c? < 2 is excluded from the outset due 
to the absence of ordering in the bulk. For 2 < d < 4 
a very important solution for the correlation function at 
criticality was obtained, however, only for the "magic" d- 
dependent values u*{d) of the coupling constant u. This 
solution yielded the anomalous dimensions ?7|| = d — 2 
and ri± = {d — 2)/2, with which all other surface criti- 
cal indices could be derived from the scaling arguments 
(see, e.g., Refs. |ll|Jl^ ]). These critical indices depend 
on d and differ from the mean- field ones, as well as from 
those of the global-constraint semi-infinite SM. In the so- 
lution of Bray and Moore the required coupling constant 
u*{d) vanishes for d = 4, i.e., the model simplifies to the 
exactly solvable and trivial Gaussian model showing the 
mean-field behavior. On one hand, it seems reasonable, 
because the critical indices indeed simplify to their mean- 
field values for d > 4. On the other hand, one could 
desire more detailed information about the critical be- 
havior of the models with u ^ u* , in particular for u 
and d > 4. Unfortunately, no extension of the analytical 
results of Bray and Moore in this direction or for the off- 
criticality case, is possible. Also a numerical solution is 
hampered for this model by insurmountable difhculties. 
Apart from the obvious impossibility of handling the in- 
homogeneities on the scale of the lattice spacing with a 
continuous field theory, it turns out that this model can- 
not be solved numerically at all because near the surface, 
where the boundary condition is set, the continuous ap- 
proximation does not apply. In fact, this is an example 
of a situation in which using a field-theoretical approach 
in statistical mechanics brings only disadvantages. By 
contrast, the ASM formulated from the beginning in its 
true form on a lattice leads to the ASM equations which 
are well defined and suitable for numerical solution, and 
also can be considered continuously far from the surface. 

The main body of the article is organized as follows. 
In Sec. m the system of equations describing the ASM in 
zero field is written down. Its bulk solution, which differs 
from the well known solution for the spherical model by 



the incorporating the uniaxial anisotropy, is studied for 
all lattice dimensions. The quantity playing the central 
role in the theory, the gap parameter Gm is related to 
the reduced energy density [/„. The continued-fraction 
formalism, which is mainly needed for the numerical so- 
lution of the ASM system of equations, is described. Sec- 
tion |l| concludes with the results for the variation of Gn 
far above and far below Tr. In Sec. [II the ASM is solved 



analytically in low (1 < d < 2) and high (4 < d) lattice 
dimensions, starting from the exact solution for the one- 
dimensional "toy" model. The energy density profiles 
and spin correlation functions are calculated analytically 
in all possible cases. In Sec. IV the most interesting 
case 2 < d < 4 is investigated. Analytical solutions are 
obtained for the isotropic model at criticality and away 
from the isotropic criticality. In Sec. ^ the semi-infinite 
ASM is numerically solved in the whole range of dimen- 
sions at T > Tc- The results for the energy density profile 
and correlation functions are presented. In Sec. VI the 



main results of the paper are summarized and compared 
with the results of other approaches. 



II. BASIC RELATIONS 

A. ASM equations 

We start from the Hamiltonian of the uniaxially 
anisotropic classical D-vector model, which, in the ab- 
sence of the magnetic field, can be written in the form 



niairriaj , (2.1) 

a=2 j 



where va.i is the normalized Z?-component vector, |mi| — 
1, and r/ < 1 is the dimensionless anisotropy factor. In 
the mean-field approximation the Curie temperature of 
this model is Tf^^ = Jo/D, where Jq is the zero Fourier 
component of the exchange interaction. It is convenient 
to use T^^^ as the energy scale and introduce the di- 
mensionless temperature variable = T/T^^^. Using 
the diagram technique for classical spin systems p9| , ^ , 
in the limit D — > oo one arrives at the closed system of 
equations ||3^ ] for the average magnetization nii = {nizi) 
and the correlation function of transverse (a > 2) spin 
components: Sij = D{maimaj) ■ This system of equa- 
tions describing the anisotropic spherical model consists 
of the magnetization equation 



the Dyson equation for the transverse CF 



(2.2) 



(2.3) 



and the kinematic relation playing the role of the spin 
constraint on a lattice site i. 
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(2.4) 



for: 



0, 



Too = 



(2.9) 



Here Sij is the site Kronecker symbol, = Jij/Jo, and 
the so-cahed gap parameter Gi is the one-site cumulant 
spin average D {m aim ai) cum /(^ renormahzed by Gaussian 
fluctuations. We will see below that Gi can be related to 
the energy density at the site i. 

The ASM system of equations is self-consistent. It is 
instructive to compare it with the MFA equations that 
can be rec over ed via the following steps, (i) Ignore corre- 
lations in (2.3), which leads to sg = 9Gi. (ii) Express Gi 



through nii with the help of_(2.4) to get Gi — {l — rnD/O. 
(iii) Insert the latter into (2.2) to obtain the closed equa- 
tion for magnetization. The form of the latter is simpli- 
fied with respect to the general-D case because of the sim- 
plification of the Langevin function in the limit D ^ oo. 
The resulting rrii is zero above 0c = 1 and nonzero below 
Oc- (iv) With Gi de term ined, which is simply Gi = \/9 
for 6 > 6c, return to (2.3) to find the improved correlation 
function. It is clear in step (iv) that the MFA is not self 
consistent, even in the simpler case above Oc- This is the 
reaso n wh y the MFA value of dc found from the CF equa- 
tion (2.3) with Gi = 1/6 for the spatially homogeneous 
isotropic low-dimcnsional magnets is nonzero, in contra- 
diction with the result of more rigorous approaches. By 
contrast, the ASM equations are free from such an in- 
consistency and they correctly describe the dimensional 
effects in isotropic and weakly-anisotropic systems. It 
should be noted that in the "Ising limit" 77 = all the 
steps above leading to the MFA equations are exact, i.e., 
the Ising model in the limit Z? — > 00 is exactly described 
by the mean-field approximation. 

For the model with the nearest-neighbor interaction 
Jij on the d-dimensional hypercubic lattice in the semi- 
infinite geometry, it is convenient to use the Fourier rep- 
resentation in d' = d — 1 translationally invariant dimen- 
sions parallel to the surface and the site re pres entation 
in the dth dimension. The Dyson equation (2^) for the 
Fourier-transformed CF cr„„' (q) then takes the form of a 
system of the second-order finite-difference equations in 
the set of layers n = 1, 2, 00, 

26nCr„n' - 0-ri+l,ri' " Cr„_i^„/ = {2d6 / r])S„n' ■ (2.5) 

bn here is given by 

6„ = 1 + d[(77G„)-i - 1] + d'{l - X'^), (2.6) 
where Aq is given by 

d' 

d' 



1 



(2.7) 



and the lattice spacing oq is set to unity. The magneti- 
zation equation (|2.2|) takes the form 



2&„TO„ - TO„+l - TO„_l = 0, 



(2.8) 



with bn = 1 + d[G„^ — 1]. Since the layer with n is 
absent, one can use 



as the free boundary conditions to (2_^) and (pTq). If 
the interaction in the boundary layer differs from that in 
the bulk (see, e.g., Ref. ||ll|) the form of the boundary 
conditions changes. The constraint equations (2.4) can 
now be written as 



d'^ q 
(2^' 



.(q) 



(2.10) 



For n 3> 1 , q <SC 1 , and k ^ 1 , where k is the inverse 
transverse correlation length defined by (|2.34|) below, the 
second-order finite-difference CF equation (|2.5| ) simplifies 
to the differential equation for the Green's function 



— - + 2dGu 
dw^ 



-2d65{n-n), (2.11) 



where n is considered as a continuous variable, q 

y^K^ + q^, and 



Gin = Gn — G <C 1 



(2.12) 



is the deviation of the gap parameter from its bulk value 
G. The magnetization equation (2. 
form with q 



takes on a similar 



and without the inhomogeneous term. 
Bef ore p roceeding, let us consider the solution of equa- 
tions ( p.5|) for the important special variation of G„ near 
the surface, in which only Gi may differ from the bulk 
value G. In this case one can solve (2^) directly with the 
result 



d6/r] 



a 



n — n I 



a 



n+n 



a - 26ii 
2hi 



a 



(2.13) 



where a = b — \^b^ — 1, b is the bulk value of b„, a nd 
bn = b — bi. One of the particular cases of ( 2.13 ) is 
bn — 0, which corresponds to the MFA or to the high - 
dimen siona l {d > 4) lattices for n ^ 1 (see Sec. Ill C ). 
Here ( 2.13| ) simplifies to 



d6/r] 



'i—n' I 



(2.14) 



Another particular case is Gi = [2d/(2d— 1)]G, i.e., 
26ii = 1. As we shall see below, a variation of G„ close 
to this one is realized for low-dimensional [d < 2) lattices 
i n the weakly anisotropic case at low temperatures. Here 
13) reduces to 



d6/r] 



(2.15) 



For 1 — r/ <C 1 and q 1 one can define 26ii = 1 — c and 
simplify ( 2.15 ) to 



d6 



-q\n — n' \ 



-q(n+n'-2) 



(2.16) 



with f{q) ~ {q — c)/{q + c). This result coul d also be 
obtained solving the differential equation ( 2.11 ) with the 
boundary condition 
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an 



= 



(2.17) 



following from (2.9). The remarkable feature of this so- 
lution is that it becomes insensitive to the exact form of 
the boundary condition if g ^ c or g 3> c. In fact, as will 
be argued below, the limiting for ms of (t„„' with both 



signs of the surface-induced term, (|2.14| ) and ( p.l5|) , are 
realized for more general variations of G„, in which G„ 
differs from the bulk value in some localized region near 
the surface, n<n*. 

Note that the quantity c above is similar to the coeffi- 
cient in the surface-energy term that is introduced in the 
phenomenological field theory of phase transitions and it 
defines the extrapolation length c = (see, e.g., Ref. 
[01). This term was used, in particular, by Bray and 
Moore Q|, who have set c = oo for the ordinary phase 
transition to remove the uncertainty. We shall see, how- 
ever, that the weakly anisotropic models with d < 4 the 
microscopic solution is characterized by effective values 
of c of order k <C 1, i.e., by large extrapolation lengths. 

The equation for the longitudinal correlation function, 
cr^^,, is not coupled to the ASM system of equations, 
since fluctuations of the (only one) longitudinal spin com- 
ponent make contributions of order 1/D to the physical 
quantities, which disappear in the spherical limit. Above 
Tc in zero field this equation has the form (2.5) with rj — 
1. The latter amounts to replacing q ^ q 



in ( 2.11 ), where is th e inverse longitudinal correlation 



length defined by ( 2.48| ) below. Thus both (t„„' and 



are given by the same function of different arguments, 
and the latter is more convenient, since its argument 
spans the wider range, starting from zero at c riticality 
(kz = 0). In this limit, q^ — 0, formula ( 2.16| ) for cr^^ 
reduces to the expression 

(T^{;,'^2dd{n-l + K), (2.18) 

where the extrapolation length Ag is given by Ae = 1 /c. 

B. Energy and susceptibilities 

For t he f erromagnetic model described by the Hamil- 
tonian (2.1), the energy corresponding to the ith site, Ui, 
is determined in the spherical limit D — > cxd by the spon- 
taneous magnetization = {rrizi) and the transverse 
CF s„: 



(2.19) 



It is convenient to consider the reduced energies U = 
U/\Uo\ where Uq — —Jq/'2. is the ground-state energy. 
In the semi-infinite geometry the reduced energy corre- 
sponding to any site in the nth layer can be written as 



(2.20) 



where Unn is due to the interaction with one of the neigh- 
bors in the same layer and Un,n±i is due to that with the 



neighbors in the adjacent layers. The terms of ( 2.2C ) can 
be represented through magnetization and the layer CF 
cr„„'{q) as 



" 2d 



d'^'q 
(2nf' 



Un,n±l — — 



2d 



mnmn±i 



f7Aqf7„„(q) 
d'^'q 



(2.21) 



7Wn,n±l(q) 



In practice, only the tot al en ergy of the site Un is needed. 
Using the CF equation ( ^.5[) to eliminate an,n±i a nd th e 
magnetization equation ( |2.8| ) to eliminate mn±i in (^-21), 
one comes to the remarkably simple result 



Un = e- 1/G„. 



(2.22) 



The deviation of the energy density from the bulk density 
U is given by 



Ui„ = Un-U^ Gm/G^ 



(2.23) 



where Gi„ is defined by ( 2.12| ). This formula provides 
additional physical interpretation of Gi„, besides that 
following from the role it plays in the correlation func- 



tions (see Sec. II C ) 



The susceptibilities of a ferromagnet are related to the 
correlation functions. In the semi-infinite geometry the 
generic susceptibility is that describing the response of 
the spin polarization in the 7ith layer to the dimensionless 
magnetic field h = JqH in the n'th one. In the region 
above Tc which is considered throughout the paper, it 
can be written as 



Xann' — ^(^an)/^^Qr; 



,i0)/9. 



(2.24) 



Here, in the anisotropic case the longitudinal {a — 1 — z) 
and transverse (a ^ 1 = z) susceptibilities are different. 
For the transverse susceptibility the corresponding layer 
correlation function (t„„' (0) (the index a is dropped for 
conve nien ce) is determined by the system of linear equa- 
tions (2^). The longitudinal layer CF a^^,{0) satisfies 
the same system of equations with rj ^ 1 i n (2.6). The 
most important of the susceptibilities ( 2.24 ) is Xaii, cor- 
responding to the boundary layer. Whereas for rj < 1 
the transverse susceptibility xu is noncritical, the longi- 
tudinal one Xzii shows critical behavior with the critical 
index 711. One can also consider the response in the nth 
layer to the homogeneous field. The appropriate suscep- 
tibilities are given by 



dha 



E 



1 



(0). (2.25) 



C. Bulk limit and continuous dimensions 



In the homogeneous case, inj = m and Gi = G are 



constants, and the equation (2.3) can be easily solved 
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with the help of the Fourier transformation. This yields 
the Fourier-transformed transverse CF 



s(k) = 



6G 



1 - TyGAk ■ 



(2.26) 



The longitudinal CF Szz{^) is given above Tc by the same 
expression with r/ = 1 . Now the autocorrelation function 
Sii can be expressed as 

s^^ = = 0GP{7iG) = \-m\ (2.27) 



where 



(2.28) 



is the lattice Green's function. The quantity Ak = JiJJp 
is given for the nearest-neighbor (nn) interaction by (2.7) 
with d' ^ d. The total wavevector k is related to q above 
by k = kzBz + q, where q • =0. The last equation in 
( 2.27 ) together with the equation 



m(l - G) = 



(2.29) 



following from ( p.2[ ) in the homogeneous case completely 
describes the ASM in zero magnetic field H. The ho- 
mogeneous ASM equations for iJ 7^ can be found in 
Refs. The lattice integral P{X) has the following 

properties: 



P{X) = { 



i + xy{2d), 

l/(l-X2)l/2^ 

(l/^)ln(8/Xi), 

Wi-cXi ln(c7Xi) 
Wd~cXi, 



X< 1 

d=l 

Xi < 1, d^2 

Xi < 1, d = 3 

Xi < 1, d^A 

Xi < 1, d > 4, 



(2.30) 



where = 1 - X and W3 = 1.51639, W4 = 1.23947, 
= 1.15631, etc., are the Watson integrals. Since 
above the constraint equation (2.4) yields su = 1, 
equation ( ^.27 ) determines the value of G which increases 
with decreasing temperature. The high-temperature 
asymptote of G is G ^ e-^[l-Tfe-^ /{2d)], 6I > 1. This 
results in [/ = —if_[{2d6) for the energy in the bulk, 
which is given by ( p.22| ) with n — > 00. The criticality is 
determined by G = 1, which corresponds to closing the 
gap i n the longitudinal correlation function S2z(k) [see 
( 2.26 )]. This is the reason for calling G the "gap pa- 
rameter". Below 9c one obtains G = 1 from (2.29) and 
then 



Oc = i/P{v) 



(2.31) 



from ( 2.27 ). Here the value of the Curie temperature 
9c generalizes the well-known result of the spherical 
model 9c — l/W |25| for the anisotropic case. 



The influence of the anisotropy on the ordering in the 
ASM is rather essential. The anisotropic gap in the 
transverse CF s(k) prevents long- wavelength excitations 
(transverse fluctuations) from destroying the lo ng-ran ge 
order in two dimensions, and 9c determined by (2.31) is 
flnite for rj < 1. Moreover, the phase transition at flnite 
temperature occurs even in the one- dimensional ASM. 
This surprising result is due to the switching off of the 
longitudinal fluctuations in the limit Z? — s- 00, which are 
responsible for the breakdown of the long-range order in 
one dimension. 

Th e bulk solution of the linear system of equations 
(^.5D has the form 



<^nn' (q) 



d9 



1, 



(2.32) 



where b is given by ( |2.6| ) with G„ G. This result could 
also be obtained by the integ ration of the bulk trans- 
verse CF s(k) given by ( ^.26| ) over k^- For the weakly 
anisotropic ASM, 1 — <C 1, at small wave vectors the 



transverse correlation functions in the bulk, (2.26) and 
(^.32 ), have the form 



s(k) 



2d9 



(q) 



d9 



where k is defined by 



= 2d[l/(ryG) - 1] =^ 2d[l - 77G] < 1. 



(2.33) 



(2.34) 



One can see that the transverse correlation length = 
1/k increases without diverging with decreasing temper- 
ature down to 9c and remains constant below 9c, in ac- 
cordance with the behavior of G. 

The field-theoretical multiple-component 0'* model 
used by Bray and Moore extends in a natural way 
for arbitrary noninteger lattice dimensions d. The dis- 
crete structure of the lattice which is important near the 
surface is, however, lost in such a model. A better way 
to get a continuous-dimension model to study crossover 
between different lattice dimensions is to consider the d' 
translationally invariant dimensions as continuous, pre- 
serving the dimension z perpendicular to the surface as 
discrete. This amounts to making the long-wavelength 
approximation 



1 , 

Ak = ^cosfc^ + — Aq, 



A' 



1 



in the whole Brillouin zone for the part of the expression 
for the exchange integral Jk = JoAk- The natural hy- 
percubic cutoff \ki\ < n and the corresponding density of 
states are modified for the q components according to 



2^---^A^./. ^'^'^ ■ 



(2.36) 



with A = y/2{d + 1) . One can check that the sum rules 
J/ M 



Ak 



rdh_d/_ r 



(2.37) 
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are satisfied. Now using ( 2.32 ), one can, instead of ( 2.27 ) 
write 



4f I '^^'i" 



1 



0GP{r]G), (2.38) 



which is the definition of P{X) in our continuous- 
dimension model. The re sulting P{X) posesses the same 
general properties ( ^.30|) . The Watson integrals W for 
some values of d are W2.5 = 2.527059, W3.0 = 1.719324, 
W4.0 = 1.321825, and W5.0 = 1.192848. 

For both hypercubic and continuous-dimension lattices 
the singular behavior of the integral P(X) for k <C 1 is 
described by 



P{X) 



W-CdK'^-^, 2 < d < 4. 



(2.39) 



Here Cd — Ad 'X- dMd, where the non-universal factor Ad 
reads 



Ad 



Sd'/i'^n)'^ , hypercubic lattice 



d'/A'^ 



continuous dimensions 



(2.40) 



and Sd = 27r''/2/p(gj/2) is the surface of the d- 
dimensional unit sphere. The universal quantity Md 
which will be needed below is given by 



Md 



VTTF 



^2~d 



27ri/2 



-r 



1 



2-d 
2 

(2.41) 



for d < 2 and 



Md 



dyy' 



d'-l 



r(d-i) 



d- 3 



cos(7r^) 2'i-ir2(d/2)^ 



(2.42) 



^2-d 



in (|2.4lD is for k < 1 



for 2 < d < 4. The factor 1 
close to unity if d is not close to 2. It is needed for 
2 - d < 1 to give 



Md = 



1 



2-d 



= ln(l/K) 



(2.43) 



with logarithmic accur acy. For d ^ I one obtains Cd 
1, in accordance with ( 2.3Q| ). 

In the anisotropic case 77 < 1, the value of G deter- 
mined from the equation 9GP{riG) = 1 approaches 1 
l inear ly just above the Curie temperature 9c given by 



1 - G - r//(r,), 



-1, 



(2.44) 



where 



liX) ^ 1 



XP'(X) 

~p(xr^ 



P'W-^. (2.45) 



For the weakly anisotropic model this solution is valid in 
the narrow region defined by 1 — G ^ 1 — 77, i.e., below 
the crossover temperature t* — {\ — ri)I{ri). For different 
lattice dimensions r* reads 



1, 

l/ln[l/(l-,7)], 
(1 - r7)(''-2)/2^ 

(l-77)ln[l/(l-r7)], 

1-^, 



d<2 
d = 2 
2 < d < 4 
d = 4 
d > 4. 



(2.46) 



For T ^ T* one has 

02/(2-rf)^ d<2 

exp{-A^^/e), d = 2, 2.0 



1-G 



2/{d-2) 

r/ Inr, 
r, 



2 < d < 4 
d = 4 
d > 4, 



(2.47) 



where, according to ( 2.40 ), A^^ = and ^2 ~ v^- 



Here the result for d < 2 is valid for Q <^ 1, i.e., a weakly 
anisotropic system can be close to criticality (1 — G ^ 1) 
in a temperature range extending far above 0c ^ 1- For 



d > 2 the Curie temperature is not small, and ( 2.47 ) 
requires t ^ 1. 

The longitudinal CFs are given by the same formulas 
( p.33D and ( |2.34| ) with = 1 . The longitudinal correla- 
tion length = l/'Sz, where 



2d[l/G-l]^2d[l-G]<l, 



(2.48) 



diverges at in different way s for the isotr opic and 
anisotropic models according to ( ^.44| ) and ( 2.47 ), respec- 
tively. The critical behavior of the ASM is, for 77 < 1, 
in all respects analogous to that given by the mean-field 
approximation. This is due to the suppression of th e sin - 
gular ity o f the lattice Green's function P{X) [see ( 2.3C ) 
and (2.39)]. For 1 — 77 <C 1 far enough from i-e., for 
T ^ T*, the system behaves isotropically and, in par- 
ticular, — Cca- Crossover at t* is analogous to that 
between the Heisenberg and Ising universality classes in 
the weakly anisotropic Heisenberg model. Here one has 
the crossover between the spherical and mean-field uni- 
versality classes instead. 



D. Continued-fraction formalism 

The linear homogeneous second-order finite-difference 
equation 



2&n,-2'n 



2„-i = 0, 



(2.49) 



which corresponds to the CF equation (2.5), has two lin- 
early independent solutions, X„ and /C„. They can be 
chosen so that X„ 00 and /C„ — > for 7^ — > 00. The 
solu tion of equation (|2.5| ) with the boundary condition 
.91) can be expressed through X„ and /C„ as 



7 



2d9 261X1 - 12 
2d6 J I„/C„/, n < n' 

where the Wronskian W„ is given by 



(2.50) 



(2.51) 



It can be shown with the help of (2.49) that Wn+i = Wn, 
i.e., yVn is independent of n. It is convenient to redefine 
X„ by replacing it by its linear combination with /C„, so 
that the redefined I„ satisfies the additional requirement 
Xo = in the non-existing layer n = . This entails 
261X1 — X2 = 0, i. e., th e first term in (2.50) becomes zero. 



The solution ( 2.50 ) can be rewritten in the form of a 
continued fraction, which is appropriate in particular for 
numerical calculations. In terms of the functions and 
a'j determined by 



X71 — 1 /X^ 



ICn-i/ICn = 26„ - (2.52) 



the solution ( 2.50 ) for n ~ n' becomes 
2de 1 



77 26„ 



(2.53) 



The functions a„ and can be found from the forward 
and backward recurrence relations 



1 



1 



ttn+l 



26„ 



26„ 



(2.54) 



being the consequence of (2.49). The initial condition 
for the first one is ai = 0. Far from the surface a„ 



approaches the bulk value a of (|2.32| ). The backward 
relation for a'j starts from a far from the surface. For 
numerical calculations in the isotropic case for 2 <rf < 4, 
a refined asymtote is needed (see the end of Sec. pV A ). If 
the denominator of ( 2.53 ) becomes zero for q = (which 
is the cas e for is otropic model at criticality for 2 < d < 3; 
see Sec. IV A ), then with the help of ( 2.54 ) one can 
obtain the relation 



1 



(2.55) 



The general solution ( 2.50 ) can b e represented through 
the diagonal Green function ( 2.53| ) via the relations 



^n,n—rn — ^n — m.n — ^n — 'm-\-1^7i — rri+2 ■ 



or, alternatively. 



n+m — l^n+m~2 ■ 



(2.56) 



(2.57) 



The consequence of these two relations is the useful for- 
mula 



(2.58) 



It is convenient to introduce the deviations from the 
bulk values 



and 



air 



6ln = 6 - 6„ = 



TJ \G Gn 



VGG^ ' 



(2.59) 



(2.60) 



where a, 6, and Gin are defined by ( p.32| ), (2^), and 
( ^.12 ), respectively. The recurrence relations for the de- 
viations ain and ai„ have the form 



ai,n+i 



and 



a. 



a^(26in -t- ain) 
1 - a(26i„ -I- ain 

a^(26i„ -Hq;i„) 
1 - a(26i„ + q:'i„ 



ail 



-a 



(2.61) 



0. (2.62) 



In terms of the deviations ain, ai„, and 6i„ the CF ( 2.53| ) 
can be now written as 



de 



where 



'{ain + ai„)/2 - 6i„ 



(2.63) 



(2.64) 



plays the role of the self-energy part for the s pin CF. 
At and above 9c the constraint equation ( ^.10 ) can be 
rewritten as 



where the bulk result is given by (2.63) without E„. 

Since for n ^ 1 the quantities ain and a^of (^.59^ ar e 
small and n in the recurrence relations ( |2.6l| ) and ( |2.62| ) 
can be treated as a continuous variable, these relations 
can be reduced to the fir st-order nonlinear differential 
equations which for q = \/ + q^ <^ 1 have the form 



— ai„ = -2gai„ + a-^ 
dn 



26i„, 



-^a'ln = -29a'i„ + {a'lnf + 26i„. 



(2.66) 



These Riccati equations can be transformed to linear 
sec ond-o rder differential equations which are equivalent 
to (^Hl) . 



E. Variation of the gap parameter at low and high 
temperatures 

The ma in prob lem with the solution of the ASM equa- 
tions (2.5)-(2.10) is to find the variation of the gap pa- 
rameter Gn that plays a fundamental role in th e theory. 
Its inhoniogeneous part Gin defined by (2.12) is anal- 
ogouos, as we shall see, to the function V{z) with the 
opposite sign, which was considered by Bray and Moore 
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||48| . The inhomogeneity Gin results from the deficit of 
interacting neighbors in the region near the surface and 
is positive. The simplest case in which Gn can be found 
analytically is T = 0. Here, for the magnetization one 
has m„ — 1 everywhere, and G„ determined from ( |2.8| ) 
reads 




, n ~ 1 
n>2. 



(2.67) 



This result also shows that the boundary layer, n = 1, 
is distinguished from all other ones. This feature that 
is beyond the scope of the continuous field-theoretical 
approaches can be observed in the whole temperature 
range. In particular, at high temperatures (0 3> 1), or in 
the whole region 9 > 9c in the "quasi-Ising" limit 77 ^ 1, 
the variation of G„ can be found with the help of the 
high-temperature series expansions (HTSE) . This can be 
most conveniently done u sing th e diagram technique for 
classical spin systems [^9 50 5^,Q . The result for the 
hypercubic lattice has the form (77/6* <C 1) 



Gir 



n = 1 



i f-l-Y 

9 \2d9) ' 

4(d-l) / 7? x2(n+i) 



(2.68) 



The terms with tt. > 2 are slightly different for the 
continuous-dimension model. One can see, again, that 
the boundary layer is distinguished. For the layers with 
n > 2 the expected leading diagrams of order 2n are 
cancelled, and the result is much smaller. 

Calculation of Gin in other cases requires more spe- 
cialized methods, which will be considered below. 



III. SURFACE-INDUCED CORRELATIONS IN 
LOW AND HIGH DIMENSIONS 

A. Solution for the "toy" model d = 1 

In on e dim ension the solution of the ASM equations 
(U) ~ ( pTO| ) is greatly simplified since d' = and there 
is no integration over the wave vector q. The quantity 
hn of ( |2.6| ) reduces in this case to &„ = l/{riGn)- For 
the autocorrelator s„„ one has si mply s„„ — Unn ■ Above 
9c the constraint equation (2.10) becomes (t„„ = 1, i.e., 
all 

ann are equal to each other. This means that for the 
one-dimensional ASM the transverse susceptibilities with 
respect to the layer field, x±nn = o'nn/9 = 1/9, are the 



same for all layers. Now, from the relation (2.58) follows 
an+i = ce '„- Us ing the latter together with the recurrence 



relations ( 2.54 ), one can write the constraint equation as 



1 ^nn. 



261 f l/[2&i - (26i)"i], n = 1 
7] I l/[a,7^ - an], n>2. 



(3.1) 



This implies that all a, 
equal to 



and a'n, except for ai = 0, are 



a = [v/ry2 + 9^- 9]lij, 



(3.2) 



which is with the bulk value given by (2.32). Then for 
n> 2 with the help of ( ^.32 ) one can identify bn — b = 



[a + a-^]/2 = v/772 + 9'^/r] and a/P^ = 9/r]. The 
boundary-layer quantity bi can be determined directly 
from (|t|) with the result 2bi = = + 9^ + 9]/t]. 
Now the exact result for G„ = 1/ {rjbn) can be written in 
the form 



Gn — 




(3.3) 



with 9c = \J\ — rf- ■ One can see that Gi > G„>2 = G. 
In particular, in the weakly anisotropic case, 1 — 77 <C 1, 
at criticality one has Gi = 2(1 — k), which is nea rly tw o 
times greater than the bulk value G = 1 [cf. ( 2.67 )]. 
Variation o f G„ above belongs to the class s tudie d at the 
end of Sec. II A, and thus cr„„' is given by ( ^.13[ ). In our 
one-dimensional model, however, one has a — 26ii — 0, 
and the inhomogeneous term vanishes. Thus one arrives 
again at the result tT„„ = 1, which can serve as an inde- 
pendent check of the calculations. 

Now we consider the longitudinal CF cr^^, and the cor- 
responding susceptibilities. The solution of the finite- 
difference equation (p.5|) with rj = 1 and Gn given by 
(0)^has the form ( ^.13[ ) with a^^ 
(B) ], 6, EE 1/G, and b.u 



Using (3.3) one can write, explicitly. 



1/G - 1/Gi [cf. (I2.6OD ] 



1 fc f 



= ^l + 9^-9^c- 

z — 262II 



26. 



9+^/W~^ 



(3.4) 



In contrast to the transverse CF (t„„ given by ( 2.13 ), 
the inhomogeneous term in cr^^ does not disappear. For 
the weakly anisotropic ASM in the range 9-^1, the 
expression for cr^^ simplifies to 
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Kz 



[cf. (2.16)]. where k and Kz are given by 



K,-^9^^c- 



(3.5) 



(3.6) 



One can see that here the extrapolation length Ae — 
1/k ^ 1 is large on the scale of the lattice spacing. Well 
above 9c, where = k, there is no difference between 
the longitudinal and transverse CFs: cr^f^ = cr„„ = 1. 
Near 9c one has ^ k, and (3.5) shows the dependence 
on the distance from the surface. Whereas the bulk CF 
(n = 00) di verg es with the exponent 7 ,^""^ = ^ (see the 



end of Sec. UC), the semi- infinite CF (3^) does not for 
any finite n. In the boundary layer, crfj^ takes on the 
exact form 



261 



'11 



9+^W^c 



(3.7) 
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in the whole range oi rj. It varies from 1 at 3> 
to 2 a.t 9 = 6c- At criticahty the longitudinal surface 
susceptibility is two times greater than the transverse 



one. One can see that 711 



as in the MFA. 



B. Low dimensions, 1 < d < 2 

For d > 1 the ASM equations become nontrivial be- 
cause of the integrati on ov er the wave vector q in the 
constraint equations (2.10). The deviation Gi„ of the 
gap parameter from its bulk value is now nonzero for all 
layers, n < 00. For 9 <^ 1, the system is close to criti- 
cahty, and the inverse transverse correlation length k of 
(2.34) is small and related to 6 by 



9^l/P' 



.2-d 



/Cd 



(3.8) 



[see (|2.27D and ( p.39D , cf. (|3j)]. The variation of G„ is 
for 1 < d < 2 close to that for the one-dimensional model 



above [see (O)] and can be searched in the form 



Gn 




(3.9) 



The correction terms Gi„ will be shown below to be pro- 
portional to {d — 1)9. [ Note that the definition of Gn 
here differs from that of ( ^.12| ).] 

For the variation of the gap parameter above, the spin 
CF a„n is for not too small wave vectors determined by 
the boundary condition at the surface and given to the 
zeroth order by (2.15). This makes it possible to find Gi„ 
perturbatively in the range Kn <C 1. After that one can 
study the corresponding corrections to (t„„ , which prove 
to be small for g 3> k or d close to 1. 



1 . Variation of Gn 

The quantities Gi„ for n > 2 can be found from 
the constraint equation (2.1C) written in the form 
(d'/A''') /dg/^iA2(7„„ = 0, where 



fn+l,n+l ~ 2ct„„ + Cr„-l^ri-l- 



(3.10) 



The s econd-order difference above can, with the help of 
(2.58), be rewritten as 



AV„„ = 2[-2a5i„ + {y/b^ - 1 - bin) 

X {ain + a'ln) - Oiina[n]<^nn- 



(3.11) 



Here the first term is the only one that is important in the 
long-wavelength region, q ^ k. By integrating it one can 
set a 1, which yields simply — 46i„ = — 4dGi„ in the 
constraint equation. All other terms make a contribution 
from (7 3> K, and the only important one among them is 
the term containing (6^ — l)^/^ai„ where ain is induced 
by the surface. The latter can be found by comparing 



(^.15|) with ( 2.6S ), where the small terms a[j^ and bin 
are neglected. This results in 



air. 



2(62 _ l)l/2^2„-l 



2q 



1 + 



exp{2qn) + 1 ' 



(3.12) 



where the second form is valid for g ^ 1. Using ( 2.15| ) 
for q^ K, one finally obtains 



Gm^l^ j\qq''-\l~a^)a^i-'\ 



(3.13) 



where the relation 2a\/b'^ — 1 = 1 — has been em- 
ployed. This expression is explicitly small for 9 ^ k^^'^ <C 
1 or for d' = d — 1 < C 1, as was said above. For n ~ 1, 
integration in ( |3.13 ) extends over the whole Brillouin 
zone, and Gi„ is nonuniversal. It decreases with n since 
a < 1. In the range 1 <C n <C 1/k the integration 
is cut at q ^ '^/n, and one can use y/h^ — I = q and 



v2(n-l) 



-2qn 



. Expression (3.13) then simplifies to 



Gin 



d' 9T{d) r(d) «;2 



A''' {2n)'^ dMd {2n)d' 



(3.14) 



where the second form is explicitly unive rsal. Here 
is given by ( ^.41 ), or, near d = 2, by ( 2.43 ). In two 
dimensions the result above regularizes to 



Gin — 



8n2 ln[l/(aK)]^ 



(3.15) 



where a '--^ 1 is a nonuniversal factor. For d > 2 the 
values of Gi„ are no longer small, and the meth od u sed 
above fails. At distances Kn>l the integral in ( 3.13 ) is 
dominated by q < k, where cr„„ and ai„ no longer have 
the forms ( 2.15 ) and ( ^.12 ). One can write Gi„ in the 
whole range of nn in the form 



Gin — 



r(d) K^g(Kn) 
dMd (2Kn)^ ' 



where g^nn) is a crossover function, 
asymptote of Gi„ for ktt, ^ 1 is 

Gir. 



(3.16) 



The expected 



(3.17) 



with some exponent C,. The analytical calculation of the 
prefactor here seems to be very difficult. 

Th e val ue of Gn is fixed by the constraint equa- 
tion (2.10) in the first layer. We will see below that 



Gin ~ K <C , thus Gn is simply given by 



Gn 



n=2 



9 d' 
2 A^ 



dqq'^ 



(3.18) 



where ( 3.13| ) has been used. Although the derivation 
above becomes invalid for d close to 1, the resulting ex- 
pression has a well defined form Gn = (2 in this 
region, which will be confirmed below by another method 
and is in accord with the result Gi = 2 (1 — k) in one di- 
mension [see the discussion below (p^)]. 
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2. Correlation functions near the surface 

In low dimensions the quantities Gi„ are small, and 
one can try to find the corresponding corrections to the 
correlation functions perturbatively. It is clear, however, 
that since Gi„ can be respo nsible for the gap in the cor- 
relation function [see ( 2.16| )] the direct perturbative ap- 
proach to f7„„ can be inefficient. It is more convenient to 
use perturbation theory with respect to the self-energy 
part S„ in (2.63). For the variation of G„ of the form 
(pM in the expression for the latter one has 



foil 

bin 



' dGin, 



n > 2, 



(3.19) 



and the quantities ct u, an d a[j^ can be found from the 
recurrence relations (2.61) and (2.62). These equations 



simplify, since ai„ for n > 2, as well as all «!„, are 
induced by Gin and small. For the first layer, taking into 
account an — —a one obtains in the long-wavelength 
region 



2d9 



o-ii = - 



where 



(7 + Ai((?,k)' 



1=1 



(3.20) 



(3.21) 



and the dependence of Ai on k is due to that of Gi„. The 
longitudinal CF aff is given by the same expressions with 
Qz = \/ n'i + inst ead oi q= ^ 



q ■ 



In the sum (3.21) only the term with Gn is unknown 



and its value follows from th e con strai nt equ ation (2.65) 
which can, with the help of ( 2.38| ) and (2.39), be put into 
the form 



/ dqq''-' 



q + Ai{q,K) q 



0. 



(3.22) 



For 1 < d < 2 the integral here is dominated by q ^ k 
and its upper limit can be set to oo. This implies that 
Ai K, i.e., the individual terms of the sum (3.21), each 
being of order 6 ^ k^^"^ ^ k, are nearl y can celled. The 
ensuing expression for Gn is given by ( 3.18| ) above. 

As a result of the cancellation of the leading terms in 
( 3.21 ), taking the next terms into account may become 
necessary. The mth-order terms in Gi„ are proportional 
to and they are small in comparison to k <C 1 for 

d < 2 — 1/m. In particular, the first-order perturbation 
theory written above neglects the terms starting from 
m ~ 2 and it is, in general, valid only for d < 1.5. The 
second-order perturbation theory works for d < 1.75, etc. 
The solution of the problem seems to undergo an infinite 
number of crossovers for d approaching 2 and it should be 
rather complicated. Analytical solutions can be obtained 
for d close to 1 and for g 3> k, as well as in the marginal 
case d = 2, with only logarithmic accuracy. 



For d' = d - 1 < 1 one has Md ^ l/ {d - 1), and the 
main contribution to the integral in ( 3.25 ) stems from the 
region q <^ n. For d — 1 the obvious solution of equation 
( |3.22D is Ai(k, n) = k. Since for d = 1 aU Gi„ with n>2 
disappear, this leads to the one-dimensional result Gn — 
— k/2. In the first order in d — 1 one can still neglect the q 
de pende nce of Ai(g, k). Then the perturbative solution 
of ( 3.22| ) leads to the simple result 



Ai(k, k) = dn, 



1 < 1. 



(3.23) 



On the other hand, the sum ( |3.21 ) near d = 1 consists of 
two contributions: 



Ai(.,.) = -2dGn-^E]^ 

1=1 

= -2dGii-K2-^(l-K''-i), 



(3.24) 



with logarithmic accuracy. Comparing it with ( 3.23| ) 
yields 



Gil — 



1 + (d - 1) In - 

K 



(3.25) 



For the analysis of the limit q ^ k and of the behavior 
of the longitudinal CF a^^^ for small wave vectors, it is 
convenient to represent the quantity Ai(g, k) above as 



Ai(g, At)- Ai(0,K) + Ai(g,K), 



(3.26) 



where 



Ai(0,K) = -2d5^Gi, 
1=1 

oo 

Ai(g, K) = 2dJ2il-e-^^')Gii. 



(3.27) 



1=1 



The quantity Ai(0, k) determines the gap of crfj^ at criti- 
cality {(jz — 0) and it is related to Ai(k, k) studied above 

by 

OO 

Ai(0, k) = Ai(k, k) - 2d^(l - e-2-')Gi, 



1=1 



Ai(k, k) 



KT{d) 



dx 



d Jo ^ 



-{l-e-'ngi^), (3.28) 



where the function g{Kn) was introduced in (3.16). The 
integral term in this formula stems from the region x = 
HLn 1 where ginn) is unknown. Even near d = 1, where 
Ai(k, k) is given by ( 3.23 ) to the first order in d — 1, 
calculation of this term needed to find Ai(0, k) would 
require additional efforts. In general, there is no appar- 
ent way to calculate analytically the gap in afjj . On the 
other hand, the existence of this gap can be anticipated, 
since the anisotropic model shows the mean-field critical 
behavior in all cases. At criticality for q <C 1 the longi- 
tudinal CF can, after the expansion of ( 3.27 ), be written 
in the form 
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2de 



Ai(0,k) + 



(3.29) 



where the stiffness of the longitudinal spin fluctuations 
is given by 



A=l 



nd) 



dx 



(3.30) 



which is again determined by the region ktt, ~ 1, in gen- 
eral. For d < 2 the lower limit of integration can be set 
to zero, and A is a number. For d —^ 1 the quantity 
Md dive rges, and A tends to 1. For d 2 the integral 
in ( 3.3C| ) diverges logarithmically at the lower limit , and 
this divergence compensates that of Md [see ( |2.43D ]. As 
a result, one obtains A = 2 + 0[1/ ln(l/K)]. It should 
be stressed, however, that in fact A cannot be calculated 
perturbatively, as above, if d is not close to 1. This is 
because the term of the first-order in Gi„ gives a contri- 
bution comparable to the zeroth-order result, and so do 
the terms of all orders in Gin- 

The analytically tractable case is g S> k, in which the 
sum in (3.27) is dominated by Z ~ l/q <C 1/k and one can 
use ( 3.14| ) for Gi„. Re placi ng this sum by the integral 
and combining it with ( |2.41 ) yields the result 



Ai(<7, n) - 



r(i)r(^)i-g^-'^ 



(1 - ^'^-i)g 



2-d 



(3.31) 



Here the factor 1 — q reflects the logarithmic diver- 
gence of the integral at the lower limit for d close to 2, 
its counterpart 1 — comes from Md, and the factor 
1 — k'^~^ is due to the logarithmic divergence of the in- 
tegral at the upper limit for d close to 1. For d not close 
to 1 or 2 these factors can be dropped. For d — 2 the 
expression above regularizes to 



Ai(g,K) = gln(l/<z)/ln(l/K) 



(3.32) 



with logarithmic accuracy. One can see that the q- 
dependent term calculated above is s maller than the 
leading term q in the denominator of (3.20) by a fac- 
tor [n/qY^'^ ^ 1, which justifles using the perturbation 
theory in Gi„. This term depends on k and thus signifies 
the gap in the correlation function. For 2 — d ^ 1, the 
latter has a long-tail character, which makes the pertur- 
bation scheme slowly convergent. This is in contrast to 
the bulk behavior, where this gap tail stems simply from 
the expansion oi q — \J k?- + cp- and has the fast-decaying 
form q{K/qY /2. The most drastic situation is realized 
for d = 2, where the gap tail is logarithmic and the ap- 
plicability of the perturbation theory requires fulfillment 
of the very difficult criterium \Yi{l/q) <C ln(l/K). 

One can improve the perturbation theory by taking 
into account the terms of the second order in Gi„ in the 
denominator of ( 3.20 ). These terms have the form of 
double sums over the layers and for g 3> k they make a 
contribution of order q{K/q)^^^~'^^ to Ai(q, n). For d = 2 
the formulas simplify and one obtains the contribution 



qR^ /2, where R is the ratio of logarithms in ( 3.32 ). In 
fact, in two dimensions one can sum up (with logarithmic 
accuracy) all the orders of the perturbation theory in 
Gin- This is possi ble because the equation ( 2.11 ) with 
Gin given by ( 3.15 ) can be exactly solved in terms of the 
modified Bessel functions. The corresponding calculation 
will be presented below; here we discuss some further 
features of the semi-infinite ASM for d <2- 

In each order of the perturbation theory, Ai(g, k) can 
be represented in the form ( |3.26 ), where Ai{q, k) is deter- 
mined hy I 1/q ^ 1 and is thus universal. The quan- 



tity Ai(0, k) is fixed by the constraint condition ( 3.22 ), 
where g ~ k 1, and it is thus universal, too. The 
same can be shown for all values of n, i.e., the spin CFs 
in the semi-infinite weakly anisotropic ASM are univer- 
sal in the whole half-space for d < 2- In other words, in 
this case the strong scaling is realized, which manifests 
itself in the independence of the CFs of the lattice spac- 
ing oq. Alternatively, this can be se en from the fact that 
the Green function equation (2.11) with the boundary 
condition ( 2.17 ) is applicable everywhere, because Gi„ 
are small and cr„„' is a smooth function of n in the long 
wavelength region. The nonuniversality of Gi„ in sev- 
eral boundary layers does not play a role because cr„„/ 
is sensitive only to the cumulative action of Gi„ from 
a large number of remote layers, n ~ 1/g 3> 1. For 
d = 2 there are nonuniversal logarithmic corrections to 
the strong scaling. For d > 2, as we shall see, the scaling 
is realized only in the asymptotic region n ^ 1. 

The expressions for all ct„„ can be obtained recurrently 
from (Til with the help of relation (2.58), which results 



■ Oi2 



-ail- 



(3.33) 



Here, for g ~ k all and a'n, except for ai = 0, are 
close to unity. Thus all CFs in the layers near the sur- 
face are close to each other. At large distances th e sp in 
CFs should cross over to the bulk result ( 2.33| ) [cf. (3.5)]. 
For q ^ K, this crossover occurs at Kn ~ 1 and it can- 
not be described analytically, since the solution depends 
on the values of Gi„ in this region, which have not been 
determined. For the wave vectors q 3> k the equation 
for cr„„ can be solved perturbatively in Gi„ for all dis- 



tances. This can be seen if one rewrites (2.11) in terms 
of the dimensionless variable z = qn and uses ( p^ for 
Gi„. Then the term with Gi„ in the equation becomes of 
order (K/q)^~'^/z'^ and can be considered perturbatively 
for d < 2. The solution is given ( 2.16 ) with q ^ c plus 
a correction term. Near the surface, z = gn <C 1, to 
the lo west order in z this solution can be put in to th e 
form ( 3.20 ) with A„(g, k) = Ai(g, k) given by ( 3.31 ). 
For z = gn 3> 1 one obtains 



de 
q 



r(d) 



K 

2dMd [q 



(3.34) 



which is in fact the expansion of (T„„ = dQj ^cp- — 2dGi„ 
for large g. The latter result has a simple interpretation: 
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For g S> 3> K the surface term in (t„„ is neghgible, 
and the correction to the bulk result is due to Gi„ from 
the narrow region \n" — n\ <C l/g, where Gi„ does not 
change significantly (the local correction). 

For the longitudinal CF cr^Jj, the corresponding effec- 
tive wave vector = y^tj+g^ can be smaller than k, 
especially at criticality, where Qz = <?• For ^ k the 
exponential decrease of Gi„ for Kn>l [see ( ^.17 )] comes 
into play. As a result, in the range n>n* 1/k the free 
solution (2. If) is realized again. The disappearance of 
Gin for n ^ 1/k and the ensuing free solution is also char- 
acteristic for dimensions 2 < d < 4. As we shall see be- 
low, for d > 4 the value of n* is of the order of the lattice 
spacing, and the free solution is realized in a wider range. 
In contrast to the case q ^ k considered above, here the 
sign of the surface term in cr^^, is negative. One can ar- 
gue, in general, that the form of the free solution in the 
region n* <n<l/qz is (2.16) with a coefficient /{qz/n) 
in front of the surface-induced term. The plausible as- 
sumption about the form of / is / = (gz — Ceff)/(gz-l-Coff), 
i.e., the surface term of the CF changes sign as a function 
of the wave vector. As a justification one can stress that 
from the distances n ^ n* the region n<n*, where the 
form of the CF is complicated, is seen as narrow. Thus 
one can replace this region with an effective boundary 
condition of the type ( 2.17 ) set at n '--^ n*. The quantity 
Coff can be expected to be of order k in dim ensi ons low 
enough. This is the exact result for d = 1 [see (3^)] and it 
will be shown numerically to hold for d < 4. This implies 
that the extrapolation length X^, — 1 /ccS is of the order 
of the transverse correlation length ^ca, which (although 
nondivergent) is much greater than the lattice spacing 
for the weakly anisotropic systems near criticality. The 
important implication of the negative sign of the surface 
term in a^^, at small wave vectors is the gap in cr^^, at 
any finite distances from the surface, even at criticality. 
It is clear, although difficult to prove rigorously, that the 
existence of this gap in the asymptotic region, where the 
free solution is realized, also entails the gap in tr^^, near 
the surface. 



3. Correlation functions for d — 2 

Dimension d = 2 is the marginal one between d < 2, 
where the characteristic wave vectors in the spin CFs are 
g ~ K, and d > 2, where g ~ 1/n are important. The 
solution with logarithmic accuracy for d = 2 can be ob- 
tained if these ranges are separated by many decades, i.e., 
ln(l/K) ^ Inn. In fact, the solution ( 3.15 ) for Gi„ has 
been obtained under this very restrictive condit ion. W ith 
this form of Gi„, the Green function equation (2.11) can 
be solved exactly in terms of the modified Bessel func- 
tions. This will be shown in more detail in Sec. IV A, 



^(^„)l/ln(l/«) 



(3.36) 



and the result has the form (4.2) with 



ln(l/K)_ 



(3.35) 



For 



gn ^ 1 the expression for (t„„ simplifies to 



[ cf. ([4.9^ ]. This can be represented in the form of type 
(|3.20|) with An = g[l - (gn)-i/i"(i/'^)]. For not too smaU 
g one can expand A„ in powers of ln[l/(gn)]/ ln(l/K) to 
obtain the g ap tai l of the spin CF. In this way the first- 
order result ( 3.32| ) for n = 1 and all other orders of the 
perturbation theory in Gi„ are recovered. 

At small wave vectors g, solution for (T„„ is determined 
by Gin cd Kn 1 where the latter are unknown and 
hence (t„„ does not have the form above. In fact, here 
the gap A„(0, k) in cr„„ manifests itself, and it turns out 
to be much larger than the bulk gap k. The dependence 
of A„ (0, k) ca n be o btain ed fro m the constraint equation 
(^.65|) , where (|3.36|) and ( ^.33|) are used and the integra- 
tion over q is performed between A„(0, k) and 1/n. In 
this way one comes to an interesting formula: 



A„(0,k) - K'"7n, lnn<ln(l/K). 



(3.37) 



He re the critical index ln2 results from the fact that cr„„ 
of ( 3.36| ) is about 2 times greater than in the bulk; thus 
the gap in this region should be correspondingly greater 
t o sa tisfy the constraint equation. The coefficient in 
( |3.37 ) cannot be determined in the logarithmic aprox- 
imation; it should become universal for n ^ 1. This 
method is rather rough and it cannot distinguish de- 
tween the transverse and longitudinal correlation func- 
tions. One can expect that they differ by a numerical 
factor, as was confirmed by numerical calculations, which 
have been done, how ever, in a range of k not small enough 
to confirm formula ( 3.37| ) itself. 



C. High dimensions, d > 4 

As we have seen above, for the weakly anisotropic ASM 
in low dimensions the problem can be solved analytically 
for Kn 1 due to the separation of the g ranges in the 
transverse CF cr„„, as exemplified by ( 3.11 ). The surface- 
induced term is important for g ~ 1/n 3> k, whereas the 
term induced by Gi„ dominates for g ~ k <C 1. In high 
dimensions the separation of the g ranges of both terms 
also takes place, although in a different form. Whereas 
the g range of the surface term remains g 1/n (for 
Kn ^ 1), the Gi„-induced one dominates in the range 
g ~ 1. Thus the ranges separate at distances n ^ 1, 
where the problem can be solved analytically. Henceforce 
in this subsection we will consider the close-to-criticality 
case K <C 1; otherwise, at distances n ^ 1 de viations 
from the bulk values will be too small [cf. ( |2.68 )]. 

The surface term of cr„„ has in high dimensions the 
oppos ite sign, as compared to the low-dimensional result 
( ^.15 ), and for n ^ 1 its form does not depend on the 
details of the behavior in the region close to the surface. 
As we will see below, for d > 4 the quantity Gi„ decays 
fast with n, and for the calculation of the surface term it 
can be neglected starting from some n* ^ 1. For n > n * 
one can use for cr„„' an expression of the type (2.13), 
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where 26ii is replaced by some quantity determined by 
the region n < n*. In contrast to the low-dimensional 
case, there is no reason for the substitute for 2&ii to be 
close to unit y, bec ause the variation of G„ is no longer 
close to (3^), or in ( |3.9| ) Gin are no longer small. At large 
distances and small wave vectors one can thus use (2.16) 
where the coefficient /(g) in front of the surface term is 
close to —1, since the quantity Ccg should be of order 
unity (the extrapolation length Ae = 1/ceff comparable 
to the lattice spacing). 

Since postulating (T„„' in the form ( 2.16 ) is not quite 
a rigorous procedure, let us consider another derivation 
based on the continued-fraction formalism described in 
Sec. 



[ID 



^ Fo r n ^ 1 one can employ the differen- 

t ial e quations ( 2.66| ) for the quantities ai„ and a'l^ of 
(2.59). The boundary condition for the second equa- 
tion is a'lao — 0; the quantity a'i„ is generated solely 
by Gin and it is not related to the surface term of ann- 
For the equation for ain the boundary condition cannot 
be set on the surface, since this equation is invalid for 
n ^ 1. Thus we use the boundary condition ain = ain* 
at n — n* ^ 1, w here a i,,* is determined by the exact re- 
currence formula ( 2.61 ) in the region n < n* . For n > n* 
one can neglect 26i„ in the differential equation for ai„, 
after which it can be linearized with respect to the new 
variable 1/q;i„ and solved to give 



2gai„. 



ai„. -I- (2(7 - ain') e-Kp[2q{n - n*)] ' 



(3.38) 



Since ain* is generated by the boundary condition at 
the surface and by 6i„, which are not explicitly small for 
n < n*, one has |ain-| ^ q for small wave vectors. It 
can be seen that ain* cannot be positive, otherwise ain 
turns to infinity at n determined by exp [29(71 — n*)] — 
Q^in* /(«!«• ~2g) = 1. T hus, ai„. < 0, and in the relevant 
region n'^ n*, qn ^ I ( 3.3S ) simplifies to the form 



air, 



2q 



exp(2gn) — 1 ' 



(3.39) 



which is independent of the behavior in the boundary 
region, n < n*. For qn ~ 1 o ne has ain ~ 1/?^; thus in 



the differential equation ( p.66| ) qain ^ afn This 
implies that the metod used here works if Gi„ oc bin 
decays faster than 1/n^. As we shall see shortly, this is 
the case for > 4. No w the surface term of <t„„ can be 
found from ( 2.63 ) and ( 2.64 ) with a[j^ = bin = 0, which 

results in (2.16) with f{q) close to —1. 

For g ~ 1 and n ^ 1 the quan tity ai„ of (3.39) 



as well as the surface term in ( [2.16 ), are exponentially 
small. Here the Gin-induced term becomes dominant. 
To find the values of ai„ and a'ln in this re gion , one 
can drop the small terms af^ and {a[^)'^ in ( |2.66 ), af- 
ter which the linear inhomogeneous differential equations 
can be solved. The solution at point n is induced by 
bin" from the interval of n" around n, which satisfies 
1/q ^ 1 <C n] thus on e can treat 6i„ in 



( ^.66 ) as a constant. The solution of ( ^.6(j ) has the form 
ain = a'ln — Oibin/^/b"^ — 1, and the quantity A„ of 



( 2.64 ) reads A„ ^ —bbi n/ — 1. The resulting cor- 
rection to (T„„ given by (2.63) is of the form 



5u. 



local 



de 



bbir 



(62 - 1)3/2 



db^™ 



bulk 



(3.40) 



where is given by ( 2.32| ) or by the first term of 

( ^.16 ) . This correction is due to the local deviation of 6„ 



from the bulk value b [see ( 2.6C )] and it could in fact be 
written for g ~ 1 without calculations. 

Now the valu e of 5i„ can be found from the cons traint 
equat ion ( ^.65 ), where (t„„ is the sum of ( ^l6| ) and 
( |3.40|) . One can see that 6i„ > is needed to com- 
pensate for the n egati ve surface term. The integration 
of the local term ( 3.40 ) extends for d > 4 over the whole 
Brillou in zo ne, q '-^ \^ and can be accomplished with the 
use of (2.3j). This results in 



A^'X (62-1)3/ 



^ = ^-fP{^G)I{^G), (3.41) 



where I{X) is defined by ( |2.45| ). Since for d > 4 both 
P(X) and KX) do not diverge for X ^ 1, one can set 
X = 1 for weakly anisotropic ASM near criticality. The 
integration of the surface term is cut atg^--^ l/n<C 1 
for Kn and at g ^ g* = \/ n/n (l/rt <C g* ^ k) for 
mi ^ 1. The resulting Gi„ = 6i„/d has for d > 4 and 
n ^ 1 the form 



(-fin ~ 



Ad' P(l)/(1) 



(3.42) 



where K^{x) is the Macdonald (modified Bessel) func- 
tion . F or the hypercubic lattices the first frac tion 
in ( 3.42 ) should be replaced according to ( ^.36 ) by 

One can see that the form of 



id' 



(2.40)] 



Gin is nonuniversal. The limiting forms of ( 3.42 ) are 



Glr, 



and 



d' dT{d~2) 1 
P{1)I{1) {2n)'^- 



Gin — 



d' dr(^) K 



Kd'2P{l)I{l) ' 



Kn<l, (3.43) 



Kn>l. (3.44) 



If d is close to the marginal value d — 4 the contribu- 
tion to the integral ( 3.41 ) from small wave vectors be- 
comes large, and the separation of the g ranges in (t„„ 
no longer takes the place. Nevertheless, the problem can 
be solved analytically with logarithmic accuracy atdis- 
tances <C 1. In this case one should integrate in ( 3.41 ) 
down to g ^ 1/n 3> K where the surface term in (2. If) 
becomes important. This leads to the replacement 



P(l)/(1) 



d' 16 
A^d- 4 



1 



1 



(an 



\d-4 



(3.45) 



in (|3.43 ), a being a lattice-dependent number. For Gir, 
near d = 4 one obtains the result 
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Glr, 



16n^ {an 



id-4 



which regularizes to 



Glr, 



1 



IGn? In(an) 



(3.46) 



(3.47) 



in four dimensions. 

It can be seen that the appHcabiUty condition of the 
method used here, Gi„ = o(l/n^) for n ^ 1, is satisfied 
for d > 4. In the range 2 < d < 4 an attempt to ap- 
ply (f or K = 0) the same method yields for the integral 
( |3.41 ) a value of order ~ 6i„n^~'' that stems from the 
regi on q ~ 1/n. The integral of the inhomogeneous term 
in (i-U) is, for k — 0, determined by the same range 
of q and it is proportion al to n?~'^ . Equating both con- 



tributions according to (2.65) yields Gi„ oc 6i„ cx 1/ 



with some universal coefficient. One can see that formula 



(3.46) shows such a behavior for d < 4 where the term 
containing the nonuniversal number a is small. In fact, 
it joins smoothly the exact solution (4.1) for 2 < c? < 4 
found by Bray and Moore [Q, which will be considered 
in the next section. 

The susceptibilities of the ASM in high dimensions 
show the mean-field critical behavior. In particular, at 
small wave vectors one can drop the local contribution 
(3.40) and use ( ^.16 ) for (t„„ at distances far enough from 
the surface, n>n*. For q = and nn, K, n <C 1 both 
t rans verse and longitudinal CFs [see ( 2.16 )] simplify to 
(2.18). CFs in the surface region, n<n*, should be cal- 
culated numerically. Since an can be obtained from the 
CFs cr„„ far from the surface with the help of relation 



( 3.33 ) and the quantities a„ and are nonsingular, an 
shows the same critical behavior as in the asymptotic re- 
gion n 3> 1. The latter is characterized in part icula r by 
7ii = — ^, which can be found by expanding ( 2.16| ) for 
q = up to the second order in q = k <C 1 and using 
K CX following from (Ol and <^M). 

Finally, let us transform the CF ann'{q) to the real- 
space CF ann' (p) , where p is the distance between two 
points in the direction parallel to the surface. At large 
distances, n, n' , p ^ 1 t he relevant values of q are small 
and one can use (2.16) wit h q — q <^ c, disregarding 
the local contribution ( 3.40| ). The condition for this is 
e?" ^ [qn)'^n'^~'^, which is satisfied for d > 4 and n ^ 1. 
Then one comes to the MFA result, which, at isotropic 
criticality, k = 0, has the form Q 



O-nn' (p) OC 



1 



[p2 + („_„')2](d-2)/2 
1 

[p2 + (n + n')2]('^-2)/2 



(3.48) 



with a nonuniversal factor depending on the lattice struc- 
ture. In this expression the surface-induced term with 
n + n' , which is similar to the "image" term in electro- 
static problems, modifies its asymptotes at p 3> n, n' and 
n.-^ p,n'. These are ||,|ll| 



a„„Kp)ocl/p'^-2+''ii, 



^11 



= 2 



(3.49) 



for n, n' — const and p^ n,n' and 

(T„„,(p) « l/n''-2+„^, 77^ = 1 



(3.50) 



for p, n' ~ const and n ^ p,n' . One can see that near the 
surface correlations decay faster than in the bulk (77b = 
0), especially in the direction parallel to the surface. 



IV. DIMENSIONS BETWEEN TWO AND FOUR 

A. Isotropic model at criticality 

As we have seen above, in low and high dimensions the 
correlation function (t„„ consists of two different (surface 
and local) terms, which are dominant in different ranges 
of q. This property makes possible an analytical solution 
of the problem for kti <C 1 in low dimensions and for 
n 3> 1 in high dimensions. In the range 2 < d < 4 both 
terms are the most important in the range q ^ i.e., 
they cannot be separated from each other. Fortunately, 
the problem has an exact solution for the isotropic model 
at criticality for rt ^ 1 Q , where the anticipated asymp- 
tote of Gin far from the surface can be written as [see 
the discussion after (|3.47| )] 



Glr 



P 



2dn2 



P 



d-3 



(4.1) 



where the choice of the parameter fi will be justified be- 
low. For n ^ 1 and g <C 1 one can use the second- 



order differential equation ( 2.11 ) for the transverse CF, 
in which q = q at isotropic criticality. The latter can be 
solved in terms of the modified Bessel functions: 



'rm/If^{qn)Kf^{qn'), n< 



'nn'I^{qn')K^{qn), n' < n, 



(4.2) 



where the term C{nn'Y/'^ K^{qn)K^{qn') could also be 
added. This solution looks similar to the full discrete so- 
lution (|2.50| ), but here the constant C cannot be found 
from the b oundary condition at the surface, since equa- 
tion ( 2.11 ) is only valid for n':^ 1. 

There is, however, another method of finding (t„„/ ]4^ ] 
that avoids using the boundary condition at the surface 
and yields C = 0. The consideration starts with the 
eigenvalue problem 



P 



(4.3) 



whose solution ip[a,n) = ^/an Jf^(an) satisfies 

c 

dn ij)(a, n)il^{a' , n) — S{a — a'), 
da 'ip{a, n)ip{a, n') — 5(n — n') 



(4.4) 



and thus forms an orthogonal and complete basis on the 
semi-infinite interval. Then the Green function cr„„' can 
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be expressed through its decomposition over the set of 
eigenfunctions as 



^/ = 2d0 / da 



ipia, n)ip{a, n') 



+ g2 



(4.5) 



which results in (L2) without any additional terms. 

Now the value of /i can b e found from the constraint 
condi tion in the form ( ^.65 ) with the bulk CF given by 
(P3|), i.e., 



n \ dz ; 
'o 



,d-2 



2/^(z)A'^(z) - - 



= 0. 



(4.6) 



This integral which can be found in Ref. |Q, is zero 
for all n, if /i is given by (4.1) and 2 < d < 4. There is 
another solution, [l = (d— 5)/2 for 3 < d < 4, which leads 
to negative values of Gi„ and it should be disregarded 
for the ordinary phase transition considered here. The 
asy mpt otic form of the layer autocorrelation function cr„„ 
of (O) for gn > 1 is 



de 
q 



2 1 



(4.7) 



Here the first term is the bulk CF and the second term 
in the square brackets, dGvaJ_g^, the local contribu- 
tion analogous to that in ( p.34 ). This form of cr„„ is 
responsible for the convergence of the integral in (^^) at 
z = qn ~ 1, i.e., for q ~ 1/n 1. The latter justifies 
using the long-wavelength approximation in the scaling 
region, n 1. Here the discrete lattice structure does 
not show up in the long-wavelength behavior of (t„„' and 
in the form of Gin, thus 1/n is the only scale for q. In 
the opposite limit, gn <C 1, one can use 



2 sin(7r/i) 
to express cr„„ in the form 



[I^^iz) - I^iz)] 



ddn 
A* 



1 - 



r(l - ^) ^gnN^2M 

W+7) 



(4.8) 



(4.9) 



For d > 3 one has /i > and ann does not diverge at 
q for any finite n. Conversely, for d < 3 o ne has 
/i < 0, and thus the second singular term in (4.9) is 
dominant and it causes the divergence of (t„„ at small q. 
In the marginal case d = 3, (4.£) regularizes to 



2d0n 



In — + Co 
qn 



co-ln2-7, (4.10) 



where 7 = 0.5772 ... is the Euler constant and cq 
0.1159 is rather small. 



The Fourier transformed CF (42) looks very beautiful 
in real space H : 



2d9r{d - 2) 



1 1 



-1 (d~2)/2 



+ {n — n'Y' p^ + [n + n')^ 



(4.11) 



Here, in contrast to the MFA result (3.48), the bulk term 
and the surface-induced image term are nonadditive. The 
critic al exp onent s t^h and r]± determined analogousl y to 
d^l) and (|^ are ryy = d- 2 and r]± = (d- 2)/2 Q. 

In spite of the apparent similarity of the solution pre- 
sented here and that of Bray and Moore [Q, they are 
not completely identical. The difference is that in the 
spin vector model used here the constraint |mi| = 1 
on each lattice site is obe yed, which is accounted for 
in the constraint equation (2^). Bray and Moore used 
the phenomenological (j)'^ field-theoretical model with the 
0(00) symmetry, which has no constraint on the field 
(j). Accordingly, the self-consistent determination of the 
function V{z) |Q, which is analogous to —Gin here, is 
more complicated and can be done only for the single 
d-dependent magic value of the coupling constant u. 



A peculiar feature of the differential equation ( 2.11 ) 
i s th at its solution ( [4.2[ ) is twofold: for a given JPin 
(4.1) solutions with both signs oi jj. — ±\fi\ can be real- 
ized for d > 3 and d < 3. Accordingly, the eigenvalue 



problem (4.E) has two sets of eigenfunctions that form 
two different orthogonal and complete bases. An ap- 
parent reason for such a behavior is the singularity of 
Gin at n = in the continuous aproach, which is, nat- 
urally, not present in the original discrete formulation of 
the problem. This singularity and the concomitant loss 
of the boundary condition at the surface could be cir- 
cumvented by Bray and Moore with the eigenfunction 
trick above, which looks like a miracle (see other exam- 
ples from quantum mechanics in Ref. ||5^). However, 
the fact that the same differential equation has differ- 
ent solutions, e.g., for d = 2.5 and d — 3.5, contradicts 
the common sense and, more importantly, impedes nu- 
merical solution of this equation. The latter would be 
the only possibility in situations where no general ana- 
lytical solution is available, as in the off-criticality case; 
in that case the loss of the boundary condition creates 
insurmountable difficulties. 

The key to the paradox is that in the original discrete 
formulation there is no singularity of Gi„ , and the values 
of 2dGin are different for d > 3 and d < 3, although they 
may coincide in the scaling region, n 3> 1. The (rather 
essential) difference between the CFs cr„„' for d > 3 and 
d < 3 stems entirely from the nonscaling region, n ~ 1, 
which is not amenable to the field-theoretical methods. 
This is most p ron ounced in the limits d —^ 2 and d — > 4, 
where Gi„ of ( [4.l[ ) tends to zero in the scaling region but 
f mi of (^^) remains well defined and given by 



(4.12) 



The latter is nothing more than the particular forms of 
( ^.16 ), the difference between the two expressions being 
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completely determined by the nonscahn g reg ion near the 
surface. For d —> 2, the parameter c in ( 2.16| ) disappears 
with K in the isotropic hmit, and the coefficient in front 
of the surface terms is / = 1. For d ^ 4, the parameter 
c is of order unity, and for q ^ 1 the surface term is 
negative, / = — 1. Thus the isotropic-criticahty solution 
of Bray and Moore smoothly joins the solutions obtained 
for d < 2 and d > 4. 

To close this subsection let us loo k at how the 
continued- fraction formalism of Sec. 11 D wor ks a t 
isotropic criticality. Here in cr„„ is given by (2.63), 
1 = q in the long-wavelength region and 
■-^ can be neglected for n ^ 1. The 



where \/lP 
bin = dGi 

quantities ai„ and q;']^„ can be fo und f rom the first-order 
nonlinear differential equations (2.66) with q = q. The 
latter can be reduced to the second-order linear differen- 
tial equations and solved to give 



ain^~—\n{ e-«" [J^ {qn) + CK^, {qn)] } , 

< ^^\n{e'^^^[K^{qn) + C'I^{qn)]}. (4.13) 

Here the integration constant C should be set to zero, 
since a'^^ vanishes at infinity. The constant C remains 
undefined due to the loss of the boundary condition at the 
sur face f or the equation for ain- Adopting these results 



in (|^ and using I^{z)Kl{z) - r^iz)K^(z) = -1/z 
leads to the previously obtained expression (4.2) for n' = 
n, with the same addit ional term containin g C Thus one 



should set C = in ([4.131) . Then from (^) it can be 
seen that ai„ is no nsin gular at ^ 1 and the singular 
terms in (T„„ [see (4.9)] are due solely to a[n- Another 



way o f obtaining (4.13) is to use the definitions (2.59) and 
(2.52) to identify X„ = ^Jril^{qn) and /C„ ~ y/nK^{qn). 
The limiting forms of ai„ and a'j^„ are 



n 



2qn 



1 + — 



1 

qn 



and 



"in = - X 

n 



2qn 



1 

qn 



qn < 1 
qn > 1 

qn <s:l 
qn > 1. 



(4.14) 



(4.15) 



The results above will be used in the numerical solution 
of the semi-infinite ASM at the isotropic criticality. 

B. Away from anisotropic criticality 

The transverse correlation function (t„„' behaves simi- 
larly for the isotropic model slightly above 6c and for the 
weakly anisotropic model at or slightly above 6c. In both 
cases the behavior of (T„„' is modified in comparison to 
that at isotropic criticality due to the finiteness of the 
transverse correlation length ^ca = 1/k, where k is given 



by ( 2.34 ). For < k <C I the function Gi n h as in the 



scaling region n 3> I the form generalizing (4.1) 



Gin — 



2dn^ 



-g{Kn). 



(4.16) 



For nn ^ 1 one ca n exp ect, as is the case in other di- 
mensions [see, e.g., (3.44)], ginn) oc e"^''" with some n- 
dependent prefactor. Analytical calculation of this pref- 
actor seems to be impossible. In the opposite limit g can 
be written in the form 



g{Kn) = I — ad{KnY 



Kn <C I, 



(4.17) 



with r > and ad ^ 1. There is no guess about the 
concrete form of g{Kn) in the intermediate region and, 
moreover, even if g{Kn) is known, one would not be able 
to find a gener al analytical solution for the differential 
equation ( |2.1lD . For the field-theoretical model [||] the 
question of how to generalize the choice of the coupling 
constant u for k ^ [be it the magic value u*{d) or some- 
thing else] further complicates the problem and makes it 
quite intractable. For the ASM, however, the situation 
is not so hopeless: Some features of the off-criticality be- 
havior can be studied analytically using available small 
parameters; its general properties are well described by 
the scaling, and the rest can be done numerically. 



1. Scaling form of correlation functions 

It is convenient to start the consideration with the 
longitudinal correlation function c^^/. The latter sat- 



isfies in the scaling region n ^ 1 e quation (2.11) with 
=J> K^ + q"^ and Gi„ given by ( 4.16| ). In the generic case 
of the anisotropic criticality the longitudinal correlation 
length ^cz turns to infinity, and one has = 1 /£^cz = in 
( p. 111 ). In this case, for n = n' there are only three length 
parameters entering into the longitudinal CF: n, l/q, and 
the transverse correlation length S^ca = 1/k- Thus crfj^ 
can be written at 9c in the two-parameter scaling form 



Kn, 



y 



1. (4.18) 

K 



Away from criticality one more length parameter, ^cz, 
appears, but it does not complicate the problem. As 
can be seen from ( 2.11 ), (t^^(k, Kz, q) can be represented 
in the same form with y ^ a/ (k^./k)'^ + ■ Simi- 
larly, the transverse CF ann{i^,q) is given by ( 4.18 ) with 
y \/l + y^. For the isotropic model or above the 
anisotropi c cro ssover temperature r* [see the discussion 
following ( |2.48| )] the longitudinal CF coincides with the 
transverse CF. 

It should be stressed that in the ASM the transverse 
correlation length plays the main role, whereas the longi- 
tudinal one, which does not enter into the ASM equations 
(^.5D-(2.1C), is a subordinate quantity. This feature, 
which should be to some extent shared by the weakly 
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anisotropic classical Heisenberg model, provides a con- 
trast to the usual scaling scheme using the div erging ^cz 
as the main scaling parameter (see, e.g., Refs. |li| , [12| ). 

Let us now study the limiting forms of the scaling func- 
tion The bulk limit <i>^""'(x, y) = 1/y is clearly 
realized for z = xy = qn ^ 1. The isotropic critical- 
ity hmit <^^^°"^^{x,y) — 2xI^{xy)Kf^{xy) studied above 
is achieved if y S> 1 provided that x Both of these 
conditions imply that 1/k becomes greater than other 
length scales. For x ^ 1 the quantities Gi„ become 
exponentially small, and in the long-wavelength region 



^nn0 is given by ( |2.16[) , as in high dimensions. This 



implies the scaling function 



{x^y) 



-2{x+Xc)y' 



(4.19) 



where Xe = n\e ^ 1 is the scaled extrapolation length. 
This e xpre ssion could also be written in the form of the 
type ( 2.1(j ) , which makes no difference in the relevant 
region y ^ 1. One can see that the longitudinal CF 
at criticality does not diverge for g ^ at any finite 
distances from the surface, as in the MFA. We have seen 
above that cce = 1 for d = 1 and = (in fact, Ag ~ 1) 
for d > 4. Actually, for all values of x one has 



$(a;,0) 



2{l + x), 1 
2x, d > 4, 



(4.20) 



and the curves ^{x, 0) for all d should go between these 
straight lines. For a: ^ 1 in the wide range y <^\/x one 
can write [cf. (4.9)] 



$(x,y) 



r(l-M) (^Y'Fjy) 



(4.21) 



where fi = (d — 3)/2 and the scaling function F{y) de- 
scribes the crossover from the zero q to the isotropic crit- 
icality limit at y '--^ 1 and satisfies F{oo) = 1, as well as 
F{y) ^ y^^^ for y ^ 1. The latter requirement serves 
to kill the singularity in q for k 7^ 0; as a result (T„„(k, 0) 
behaves similarly to cr„„(0, g). For 2 < d < 3 one can 
simply use 



Fiy) = ^x « 1, y)/<l>-°"''(x « 1, y) 



(4.22) 



to find F{y). The function ^{x, 0) shows a crossover from 



$(a;,0) 



-2,1) 



for X < 1 to <l>(x,0) = 2{x + Xe) for 



X > 1 [see (}4.19|)]. 

On e can see that only the second, singular, term in 
( 4.21 ) makes cL /^-dGpGndcnt contribution to (^^^ in the 
limit q ^ 0. Specifically, for 2 < d < 4 one has 



d-2d~3 



(4.23) 



which in addition shows that the susceptibility increases 
away from the su rface , as it should. In the isotropic case 
from (|3|) and ( pl7| ) it follows that k - T^/(d-2)^ thus 



-zz,smg 



"'ii with 711 = (3-d)/(d-2) IJ. This 



result means that the surface susceptibility with respect 
to the surface field diverges in the ASM only for d < 3. 
The leading terms of Xznn near the surface are given by 



n, d > 3 

nln[l/(Kn)], d = 3 
^d-2^d-3^ 2 < d < 3 

l<d<2. 



(4.24) 



For the isotropic systems in 1 < d < 2 the bulk tran- 
sition temperature is zero. With respect to the latter, 
Xznn shows the critic al be havior Xzn n ~ with 
711 = 1/(2 - d) [see ( |2.34| ) and (^.47D]. This result is 
complementary to that for 2 < d < 4 quoted above, 
and it shows similar divergence with approaching d = 2 
from the other side. In the anisotropic case this low- 
dimensional critical behavior is reali zed in the range 
T ^ T* , where r* is given by ( 2.46 ). In the vicinity 
of 6c, i.e., T <C T*, the mean-field critical behavior with 
711 = — ^ is observed. 

It should be stressed that the critical amplitudes in 
the nonscaling region near the surface, n ^ 1, cannot be 
found in the continuous approximation. Here one should 
numerically solve the ASM equations on the lattice. On 
the other hand, it can be shown that the critical indices 
remain unchanged in the nonscaling region. The CFs 
in this region can be obtained from those in the regio n 
1 <ti n <ti 1/k with the help of the formulas of type ( |3.33 ) . 
Since the quantities an and a'^ are all nonsingular, aff 
differs from the result of the continuous approximation 
extrapolated to n = 1 by a numerical factor only. 



2. The gap tail of the scaling function F{y) 



It turns out that the form (4.17) of the function g{x) for 
X <C 1 d eterm ines the asymptote of the scaling function 
F(y) of (4.21) for y » 1, and in the region x <C 1, y » 1 
everything can be calculated anal ytically. For y 1 
the solution cr„„ of equation ( 2.11 ) at a point n <^ 1/k 
stems from the interval |n — n"| ~ 1 /q <^ 1/k around n; 
thus one can use g{Kn) in the form ( 4.17 ) and calculate 
the correction to fT„„ perturbatively in ad{nnY . The 
resulting expression for the scaling function $(x,y) of 
([4.I8D has the form 



$(x,y)^2x[/^(z)if^(z)-Q5^(z)], 



(4.25) 



where z = xy — qn, the first term corresponds to 
isotropic criticality. 



g = ad(i-M^)/y''< 1, 



(4.26) 



and the function Sp(z) reads 



Kiz)! dte 

' ZQ 



2 r+2Ai 

+ Kl{z)J^ + Iliz) 



2fi 



dte-'Klit), 



(4.27) 



with zq ^ li and = [2^r(l + /^)]^^ is a factor from 
(4.8). The part of the expression above without the terms 
containing is just what one obtains from the straight- 
forward perturbative scheme using the Green function 



18 



(4.2). The additional terms with in ( 4.27 ) can be 



rewritten in the form K'f^(z 



2 r+2fj. 



C'Z, 



i.e., they can always be added to the solution and their 
amplitude should be fixed from the boundary condition 
at the surface. Since this boundary condition is lost in 
the conti nuou s approximation, the exact form of these 
terms in (4.27) has been chosen above from the require- 



ment that the term CK'^{z), which was ruled out above 
with the help of the eigenfunction trick, does not appear 
again in the resulting expression for <I>(x,y). With such 
a c hoice one can set zq = 0, because the first integral 
in (4.27) converges at the lower limit. Now one can see 



that the terms with cancel each other, if r + 2yLt > 0, 
whereas in the opposite case they do not. For ^ < (i.e., 
d < 3) the function S^(z) can be rearranged as 

S^(z) = [2 sm{nfi)/7r]'KKl{z) + S^(z), (4.28) 

where 

E,{z) = Kl{z) dte-'il{t) + il{z) dte-'Klit), 



iliz) ^ I^iz) ~ [2sm{7rfi)/7r]^Kliz) 



(4.29) 



and 



K = 



1 2 



2 sin(7r/i) 



^r(r/2 + Ai)r(r/2 - ii)T^r/2)r-\r). (4.30) 



In ( pel ) the subtraction term with is present only 



for r + 2/1 < 0; the resulting expression is valid for both 
s igns of r + 2fi. The representation of S^(z) in the form 
(4.28) for /i < is convenient because of the cancellation 
of the divergence at terms in /^(t). For > (i.e.. 



d > 3) expression (4.28) remains valid as well, although 



the subtraction makes little sense and ^fj,{z) can simply 



be written in the form ( [4.27 ) wit h zq — = 0. 

The parameters r and aa in (4.17) should be chosen 
self-consistently to satisfy the spin-constrai nt co ndition. 
Here it is convenient to subtract equations (2.65) at and 
away from isotropic criticality from each other. Thus one 
can write 



dyy 



d'-l 



\x,y)~^x,y)]^Md, (4.31) 



where Md is given by ( ^.42 ). The integral on the left- 
hand side of (4.31) is determined by z = ~ 1, i.e., 
y 1/x ^ 1, which ju stifie s the appro ximations made 
above. With the use of (4.25) and (4.26) one can rewrite 
( 4.31 ) in the form 



2ad{j- 



dz Z 



d-2-r" 



E^{z) = Md 



(4.32) 



This equation should be satisfied for all values of a; = Kn, 
thus 



= d - 2. 



(4.33) 



Then ( 4.32 ) fixes the value of ad- 
2Md 



ad 



(d-2)(4-d) 



E= dzE{z). (4.34) 



The scaling function F{y) in ( 4.21 1 ') can now be iden- 
tified taking the limit z ^ 1 in ( [4.25 ). This leads 
to F{y) ^ 1 + 2Q K sm{Tr fi) / TT in the whole interval 
2 < d < 4. The latter can with the use of ( [4.2(j ) be 
rewritten as 



F{y)^l + Ad/y'^, » 1, 



where 



Ad = 



7ri/2 tan(7r^)r(d- 5/2) 
4(d- 2)S ' 



(4.35) 



(4.36) 



A remarkable feature of (4.35) is that the tail of F{y) is, 
for d < 4, anomalously long compared to that in the bulk, 
^buik(y) = y^/y/l + y^ ^ 1 - 1/(22/2). The sign of Ad 
is determined by fi — {d ~ 3)/2, and one has A d — for 
d = 3. This is in accord with the structure of (4.21); in 
all cases ^{x,y) is smaller than <l>'^°' "^'*(a;, y), as it should. 
In the limit d 4 the integral S of (4.34) diverges at the 
upper limit and Ad regularizes to A4 = ^. For d ^ f 
the quantity K given by ( [4.3C| ) diverges, and thus one 
can neglect S(z) in ([4.2^). In this limit Ad regularizes 



to A 



5/2 



-47ri/2/r^4) 



-0.539. The same situa- 



tion takes place for d ^ 2, where one obtains A2 — 
It should be noted, however, that for d cl ose t o 2 the 
tail of F{y) becomes extremely long [see ( 4.33[ )]. The 
validity of the present approximation for F{y) requires, 
for d ^ 2, very large values of y, which can become in- 
compatible with the condition z = <C 1 needed to 
represent $(a;, y) in the form (4.21). Actually, d = 2 is a 



special case with a logarithmically decaying gap tail (see 
Sec. pIB^ ). 

The quantity given by (4.34) is positive for | < 
d < 4 and negative for 2<d<f. Atd = f one 
has ad = Q due to the divergence of K and hence S. 
The latter could raise questions about the validity of the 
perturbation theory with respect t o adi nnY for d = f 
[should the higher-order terms in ( 4.17 ) be taken into 
account?]. As we will see below, the numerical results 
are in excellent agreement with the asymptotic behavior 
F{y) = 1 - O.539/2/I/2 for d = f and 2/ > 1. 



V. NUMERICAL RESULTS 

A. Variation of G„ 

In the symmetric phase, m = 0, the ASM equations 
were solved numerically in the following way. For a given 
variation of G„ and the value of the wave vector q i n 
]6|), the transverse CF f T„„(q ) can be found from (^.5[). 



In practice, the formula (2.63) was used, where ain&nA 



a'l^ were determined from the recurrence relations ( 2.61 ) 
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FIG. 1. Gin at isotropic criticality for different hypercubic 
and continuous-dimension lattices: (a) - general view; (b) - 
surface region, deviations from scaling. 



and ( 2.62| ). The res ult for (T„„(q) can be put into the 



constraint equation (2.4) to obtain, after the integration 
over q, the system of nonhnear equations for G„. Again, 
it is more convenient to work with the deviations from 
the bu lk va lues and to use the constraint equation in the 
form (2.65) where the subtraction is done analytically 
to avoid the loss of accuracy. The integrals over q have 
been performed in all cases over the whole Brillouin zone, 
even in low dime nsions. For the continuous-dimension 
model Isee ( p6| )] the range < q < A was divided 
into three or four log-spaced intervals (e.g., [0, 10~'*A], 
[10~^A, 10"'^ A], etc.) and the Gaussian quadratures over 
10 or 20 points were used in each of these intervals. For 
the hypercubic lattice the products of Gaussian quadra- 
tures were used, and the intergation was done with re- 
spect to the nonlinearly scaled q components Qi — q^ 
with p = 3 to redistribute the contribution of the singu- 
larity at q = more uniformly over the Brillouin zone. 
In fact, similar nonlinear transformations were also done 
for the continuous-dimension model. The resulting sys- 
tem of equations for the deviations G\n was solved with 
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FIG. 2. Gin at isotropic criticality ys lattice dimensionality 
d. The asymptotic scaling result for 2 < d < 4 and the 

asymptotic formula (3.43) are shown by dashed lines. 



a nonlinear equation solver based on Newton method. 

For the numerical solution in the semi-in finite geome- 
try, the boundary condition at n — oo in (2.62) should 
be replaced by one at some rimax S> 1. For the isotropic 
model at criticality one cannot just set a'ln^^^^ = 0, since 
a'ln slowly decays with n [see (115)]. This would spoil 
the behavior of correlation functions at small wave vec- 
tors and lead to an unphysical gap for 2 < d < 3. For- 
tunately, the asymptotic beh avior of a[j^ in the scaling 
region, n ^ 1, is given by ( 4.13 ) with c' = and it 
can be used as the boundary condition at infinity. The 
purpose of numerical calculations at isotropic criticality 
was to check the scaling solution (4.1) and to study the 
nonscaling effects at n ~ 1. The quantity fi was deter- 
mined self-consistently as a function of all Gin using the 
asympto tic fo rm of ain at n 3> 1 and q — [see the first 
limit of ( 4.14 )]. In this way the value /i = (d — 3)/2 has 
been confirmed. 

Above 9c or at the anisotropic criticality (k > 0) the 
inhomogeneities decay as exp(— 2Kn) [see (4.16)], and 
one can use the boundary condition ain^^^, = for 
«max ^ 1/k. Here the value of k should be taken rather 
sm all to study the details of the scaling function ^{x, y) 
in (4.18). Indeed, to reproduce the limit a: ^ 1 one 
should have nn ^ 1, where n* ~ 10 is the smallest value 
of n for which the continuous scaling solution holds. This 
implies, in turn, large values of rimax- In practice, calcula- 
tions were done for l — rj down to 10~®, which corresponds 
to K = ■\/2(i(l — 77) « 2.5 X 10~^ at the anisotropic criti- 
cality for d = 3. For such k the value rimax = 10000 was 
used, which corresponds to 2Kn^s_x ~ 5. Naturally, in 
this case the system of equations for Gin was not solved 
on each of 10000 layers. Instead, for n>10 only the 
"representative" layers with an exponentially increasing 
spacing between them were chosen to solve the equations. 
The values of Gi„ between them were interpolated with 
the help of the formula Gi„ = (a/n'') exp(— 2Krj) with 
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the values of a and b determined from Gin at the ends 
of the interpolation intervals. In all cases the number 
of unknowns did not exceed 50. Computations could be 
performed on a 486DX 66-MHz laptor 



The results for Gi„, as defined by ( p.l2 ) for all values 
of rt, are shown at isotropic criticality for different hyper- 
cubic and continuous-dimension lattices in Fig. 0. The 
gene ral view, Fig. 0(a), shows that the analytical result 
(4.1) is well obeyed in three dimensions outside the sur- 
face region. This result is universal and independent of 
the lattice structure; it is the same for the simple-cubic 
lattice {d = 3) and the 3d continuous-dimension lattice 
(d = 3.0). Formula (IT) has also been confirmed for 
other values of d around d = 3; the results do not dif- 
fer much from each other in the log scale and thus they 
have not been shown. I n four dimensions the results can 
be fitted with formula (3.47) with a = e'^ « 20, which 
implies significant corrections to the logarithmic approx- 
imation. In fact, the nonuniversal number a is slightly 
larger for the hypercubic lattice (d = 3), which can be 
seen in Fig. [l|(a ). In dimensions higher t han four, Gi„ 
follows formula ( |3.43| ). The coefficient in ( 3.43 ) depends 
on the lattice structure and is clearly different for d = 5 
and d = 5.0. 

Deviations from the asymptotic solution ([l.l|) in the 
region near the surface are shown in Fig. hj(b). There 
is a clear difference between the values of Gi„ for both 
three-dimensional lattices. 

The dependence of Gi„ at isotropic criticality near the 
surface on d is shown in Fig. ||. In the limit d ^ 2 the 
value Gil tends to -j^, which means that th e lim iting value 



of Gi is f , as given by the first term of (3.9), where, at 



criticality, G = 1. On the other hand, all Gi n wit h n > 2 
vanish in the limit d — > 2, in accord with ( [3.15 ), which 
disappears for k ^ 0. 

The algorithm for solving the system of nonlinear equa- 
tions for Gin based on the Newton method, which was 



used here, shows instability for d<2.3 and rimax ^ 10. 
This instability is responsible for the lack of points in 
the left part of the n = 10 curve. The reas on is that 
the integral in the constraint equation ( 2.65| ) becomes 
more and more sensitive to the region of small q where 
the integrand may become infinite due to a negative gap 
arising for some sets of Gi„ in the course of iterations. 
Away from isotropic criticality, the gap in the spin CFs 
stabilizes the algorithm. For each dimension d there is 
a minimal value of the anisotropy 1 — 77, for which the 
system of nonlinear equations for Gi„ does not show in- 
stability for n,„ax large enough, if the starting variation 
of Gi„ is chosen sufficiently close to the actual one. The 
latter is very important and necessitates using small vari- 
ations of the parameters, such as d, 1 — ry, nmaxj etc., in 
low dimensions. The minimal values of the anisotropy 
1 - ?7 are about 3 x 10"^ for d = 2.0, lO"'^ for d = 1.75, 
and 5 x 10^'' for d = 1.5. 

Contrary to the im plication of the scaling solution of 
Bray and Moore, ( |4.l| ), Gi„ do not go to zero and do not 
even show any singularity at d = 4 for any finite n. This 
is due to the correction-to-scaling terms [see Fig. |l|(b)] 
which become more and more pronounced as d deviates 
from 3. The crossover from the solution for Gi„ in th e 
range 2 < d < 3 and that for d > 4 is described by ( 3.46 ). 

Additionally, in four surface layers for the two- 
dimensional model at criticality, the dependences of Gi„ 
on the anisotropy parameter r] are shown in Fig. ^ Cal- 
culations down to l — rj — 10"'^"'^ were possible here, since 
the value of n„iax was chosen to be about 50, which 
is significantly smaller than the required Umax S> 1/k. 
The latter introduces a significant gap in the spin CFs, 
which is the artifact of cutting the ASM equations at 
"-max ^ This gap stabilizes the solution of the ASM 
equations. On the other hand, the values of Gi^ in sev- 
eral layers near the surface are pretty robust and insensi- 
tive to this defect of cr„„ . One can see that in the isotropic 
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limit, r] ~> 1, the value of Gn tends to ^, whereas all 
other Gin tend to zero logarithmically in accordance with 
In the opposite limit, 77 <C 1, the HTSE results 



(3.15 



( 2.68 ) are recovered. 

In Fig. ^ the calculated values of Gin for continuous- 
dimension lattices away from the isotropic criticality 
are represented in the scale d for m. The results for 
d = 5.0 show crossover from ( 3.43 ) to ( 3.44 ) as x = nn 



increases. A similar crossover from Gin 
Gin ~ e-2'""/n3/2 takes place for d = 4.0. The scal- 
ing is, however, not perfect here because of the logarith- 
mic corrections. For d — 3.0 the result crosses over to 
Gin ~ e~^'^^/n^-^ for nn 1. Note that there is no 
analytical solution for Gi„ in this region. 

The calculated values of Gi„ at the anisotropic criti- 
cality (or above criticality) in low dimensions (d = 1.5, 
1.75, and 2) are sho wn in Fig. [s . On e can see that the 
theoretical formulas (3.14) and (3.15) are obeyed start- 
ing from n > 10, although corrections related to the finite 
value of K and described by the function g{nn) in (3. If) 
are quite pronounced. For d = 2.0 the function g{nn) 
is nonmonotonic: apart from the exponential decrease at 
Kn 3> 1 it shows a singular positive deviation from unity 
for nn ^ 1. The fit in Fig. ^ suggests g{Kn) = 1 + 2-^/Kn 
for Acn <C 1. One can see from Fig. || that the correction 
term in g{Kn) has the negative sign for d = 1.5. The 
case d = 1.75 seems to be marginal. The values of Gi„ 
in Fig. |5| nicely follow th e dep endence Gi„ = G^^'^"^, 
where G^^^''™ is given by ( ^.14] ) with the additional fac- 
tor exp(— 2Kn) in the whole range of n. This is, how- 
ever, not an exact solution to the problem. The plot of 
Gin/Gi^'^°^ in Fig. |^ shows that for > 1 this func- 
tion begins to increase with oscillations. These oscilla- 
tions are not an artifact of cutting the ASM equations 
at the maximal layer number nmax- Numerical calcula- 
tions with different values of nmax give the same results. 
Although the ratio Gin/Gl^'^°^ is not exactly 1, its prox- 



imity to 1 in a wide range of n is remarkable, taking into 
account the strong dependence of Gi„ on n. 



B. Correlation functions 



After Gin has been determined, the spin C Fs ca n be 
f ound from (2.63) and the recurrence relations (2.61) and 
1/n^ to ( ^.62 ). The results at isotropic criticality in the scaling 



region n ^ 1, whi ch illustrate the analytical solution of 
Bray and Moore (O) for 2 < d < 4, are shown in Fig. 
|. One can see that for 2 < rf < 4 the solution satisfies 
CTrm < CTnra"^ fo^' Small wavc vcctors and cr„„ > crjj""^ for 
large wave vectors. For d = 4 one has cr„„ < (Jnn^ every- 



where, which contradicts the constraint equation (2.65). 
In fact, for d > 4 the scaling solution of Bray and Moore 
breaks down, and one has to take into account the posi- 
tive local contribution to (t„„ at g ^ 1 [see ( [3.40 )], which 
balances the constraint equation. For d = 2 one has 
everywhere, and the constraint relation is vi- 



bulk 



olated again. In fact, for d < 2 the form of tT„„ is changed 
by the gap in the region of small q, where ann < o'jjj^"', 
thus ensuring the constraint (see Fig. [lC| b elow). 

Deviations from the scaling solution ( p^.2|) in the region 
near the surface, n ^ 1, are shown in Fig. The correla- 
tion functions in different layers are related to each other 
by ( |3.33 ), where the quantities a„ and a'n are constants 
in the limit q and they approach the bulk value a 
of (|3|) with increasing n. For n 3> 1 small deviations 
of a„ and a^^ from a are responsible for the scaling form 
of cr„„ showing only a small change when n changes by 
one. By contrast, in the nonscaling region, n ~ 1, the 
spin CFs change significantly from one layer to another 
and they acquire in the range q ^ 1 nonuniversal numer- 
ical factors, relative to the extrapolated scaling solution. 
These factors, which are shown in Fig. ^, tend to 1 as 
some negative powers of n far from the surface. The ac- 
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curacy of the calculations is, however, not high enough to 
determine these powers precisely. One can see that the 
deviations from scaling are quite large an d slow ly decay- 
ing for the dimensions well above 3 [cf. ( 3.46 )]. On the 
other hand, for d well below 3 the deviations from scaling 
are mainly localized near the surface. Note the difference 
between the results for the simple-cubic {d = 3) lattice 
and the three-dimensional continuous-dimension lattice 
(d = 3.0). For the latter the deviations from scaling are 
anomalously small, which suggests the existence of an ex- 
act solution. If one assum es th at the scali ng form of tT„„ 
holds for all n, then from (3.33) and ( |2.55 ) it follows that 
a„ — \/l — l/n and a', + 1/n for q = 0. Then, 

using the divergence of ( 2.53 ) for g ^ one obtains 



6 



4+ ^l + l/n+ ^1-1/r 



3.0, (5.1) 



which is indeed a rather good approximation. It has the 
proper behavior G„ = 1 -|- l/(24n^) for n ^ 1, and 
the value Gi — 1.1082 is very close to 1.1067 follow- 
ing from numerical calculations. More careful analysis 
shows, however, that the formula above is not an exact 
solution for the ASM equations, where discrepancies of 
the type V2ln2 « 0.980 ^ 1 arise. 

At anisotropic criticality, which is the only type of crit- 
icality in low dimensions, the generic CF is cr^^. The 
transverse CF (j„„ , as well as crf^^ itself above criticality, 
can be obtained from the latter by the simple change of 
the wave- vector argument. The numerical results for tr^^ 
in the surface layer (n = 1) at anisotropic criticality in 
low dimensions are shown in Fig. |o[ One can see the 
gap a nd t he linear q dependence at small q, in accord 
with (3.29|). The values of the gap Ai(0, k) and the stiff- 



ness A in ( |3.29D determined from the fits of the numerical 
data exceed those calculated from ( 3.28| ) and ( 3.3C ). The 
reason for this is that the first-order perturbation theory 
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in Gin leading to (3.28) and (3.30) is valid for small wave 
vectors only i n the dimension range d < 1.5, as was ex- 
plained after (3.22). On the other h and, t he asymptote 
Ai((7, k) for q ^ K. which is given by ( 3.31 ), works nicely 
for d = 1.5 and d — 1.75. 

It can be seen that for d — 1.5 the results for the 
models with different anisotropies scale with each other. 
This confirms the concept of the strong scaling (down 
to the surface la yer) in low dimensions, which has been 
sug gested in Sec. IIIB 2. For d = 1.75 the scaling in Fig. 
|l0| is incomplete, which can be explained by the values of 
K not being small enough. For d — 2.0 the results do not 
scale, although either do not strongly deviate from each 
other, because d = 2 is the marginal dimension between 
the strong scaling and the asymptotic (n ^ 1) scaling. 
This behavior is illustrated in greater detail in Fig. [ll| , 
where the gap Ai(0, k) of aff is plotted as function of 
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the anisotropy. Whereas the dependen ces for d = 1.5 and 
d — 1.75 saturate in the limit of k = ■\/2(i(l — rj) going to 
zero, which confirms Ai(0, k) cx k for weakly anisotropic 
models, the almost perfectly straight line over several 
decades for d = 2.0 suggests Ai(0,k) oc Kln(l/K) in two 
dimensions. 

The numerical result for d = 2 above is quite plausible, 
because logarithms usually arise in marginal dimensions. 
This would im ply so mething like Xznn oc ln[l/(Kn)] 
for d = 2 in ( 4.24 ), where the corresponding position 
has been left empty. It seems to be, however, the third 
occasion in this work when numerical calculations sug- 
gest some qualitative features that does not follow from 
analytical considerations. For d ^ 2 the calculation 
with logarithmic accuracy leads to a nother dependence 
of Ai ( 0, k) which is given by ( 3.37 ). The applicabil- 
ity of ( 3.37 ) requires, however, such small values of the 
anisotropy that numerical calculations cannot be per- 
formed, and for larger anisotropics no other possible an- 
alytical approximations are seen. 

To shed some light on this puzzle, it is convenient to 
plot q'^~^<Jnn as a function of logg over the whole Bril- 
louin zone. The area under the curve is proportional to 
the integral over q in the constraint equation ( 2.65| ) , and 
the regions of q making contributions into the constraint 
can be well identified. Such plots show that the integral 
is dominated by q ^ 1/n for d > 2 and by 5 ~ k for d < 2. 
In the marginal case d — 2 the results for the lowest man- 
ageable anisotropy, 1 — ?7 = 3xlO~®, are shown in Fig. 
|l2| . One can see that the area under the bulk solution 
(the solid line) coincides with that under the numerical 
solutions for n — 1 and n — 20 (open circles). The 
curve for n ~ 2Q merges with the bulk curve for q > 1/n. 
Although the distance between g ~ 1 or g ~ 1/n, on 
the right-hand side, and q ^ k, on the left-hand side, 
is several decades, it is not large enough to apply the 
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logarithmic approximation, i.e., to integrate the solution 
obtained for k ^ g ^ 1/n between these limits. Also, 
K is not small enough to replac e a ( dotted line in Fig. 
0]by 1 making integration in ( 3.10| ) to obtain Gin of 
( |3.15D . The applicability condition for the formula (3.37) 
is clearly not sa tisfie d. Nevertheless, as we have seen in 
Fig. H, formula ( 3.15 ) is in reasonable agreement with the 
numerical res ults for Gin for small nn. Thus one can use 
the solution ( O ) for cr„„ in te rms o f the modified Bessel 
functions of index fi given by (3.35). This solution with 
g ^ g -I- A„(0, k) is plotted with dashed hues in Fig. |l^ , 
where the gap values Ai(0, k) = 5.2k (taken from Fig. 
|ll|) and A2o(0, k) — 2.7k were used as fitting parameters. 
The agreement with numerical results is rather good. On 
the other hand, 1/ \n {l/K) is not small enough to use the 
simplified form ( [3. 36 ) for cr„„. The corresponding curves 
deviate strongly from the numerical solution in Fig. [l^ , 
thus the final result ( 3.37 ) is not realized. And analyti- 



cally calculating the constraint integral with the gapped 
Bessel functions to obtain the simple empirical formula 
Ai(0, k) oc Kln(l/K) seems to be impossible. 

Now let us consider the numerical results for th e lon - 
gitudinal CF af^^ in the scaling representation (4.18). 
The scaling function ^{x,y) for d — 2.0 and 1.5 is rep- 
resented in Fig. 13, One can see that the asymptotic 
scaling (for n ^ 1) is well obeyed. For d = 2 in the 
surface region n ~ 1 there are small, seemingly logarith- 
mic deviations from the strong scaling, as was suggested 
above for the two-dimensional model. Here, for small x 
the results can be fitted with power-law functions, the 
exponent slowly decreasing with k. In particular, for 
1 — 1] = 3 X 10^* this exponent is 0.195, whi ch ro ughly 
agrees with l/ln(l/K) = 0.125 following from (^!36|) . For 
large and small values of y the numerical results contain 
the features of the fre e solut ion ^{x, y) — {1± e~'^^'^)/y, 



as was argued in Sec. IIIB 2. In fact, c? = 2 is a marginal 
dimension, and for d < 2 the free solution is reproduced 
for large and small y much better, as can be seen from 
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the plot for d— 1.5. The latter also confirms the strong 
scaling in low dimensions. 

The results for $(a;, 0) in the whole range of lattice 
dimensionalitiess are shown in Fig. |lj. All the curves are 
bound ed by the ones representing the exact expressi ons o f 
( 4.20 ). For small x the results are in accord with ( 4.24 ). 
The asymptotic form of the curves in the r egion a; ^ 1 , 
which is given by the limit y — *■ of (4.1E), determines 
the extrapolation length Ag. The latter is represented in 
Fig. |l5| as a function of d. One can see that for d < 
4 the extrapolation length is of the order of transverse 
correlation length ^ca, which is a "mesoscopic" length 
scale between the lattice spacing ao and the diverging 
longitudinal correlation length ^cz- 

For d > 2 the scaling form of the correlation func- 
tion holds in the asymptotic region n ^ I. The general 
view of the scaling function ^{x,y) is shown for d = 3.0 
in Fig. 16. One can figure out how the results look 
like for other dimensionalities with the help of Figs. |lj 
and lisl. The results for t he wave-vector scaling func- 

, for 2 



t ion F {y) defined by ( |4.21 
(Ol) 



< d < 3, also by 
2.5 in Fig. |l^. The dimension 



are shown for d ■ 
d = 2.5 is especially convenient since here the coefficient 
Ad in the gap tail of F{y) given by ( 4.35 ) simplifies to 



A5/2 = -47ri/Vr^(l/4) ~ -0.539. In order to reduce 
the value oi z ^ xy ~ qn, which should be small accord- 
ing to the definition of F{y), the solution for the CF in 
the first two lay ers ha s been used. Calculation of F{y) 
with the help of (4.22) yields the curves of solid triangles 
m Fig. 0. These curves do not scale with each other, 
because strong scaling does not hold for d > 2. Nev- 
ertheless, correlation functions in the asymptotic region 
n 3> 1 differ from those in the nonscaling region near the 
surface only by numerical factors, which are represented 
in Fig. ||. Inserting these factors into F{y) makes the 
results for n = 1 and n = 2 scale. These corrected re- 
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FIG. 13. Scaling representation of a^'^ for d — 2 and 
d = 1.5. Dashed lines: free solutions ${x,y) = (1 ± e^'^^^)/y 
for large and small values of y. 



suits are in excellent accord with the asymptotic formula 
d^H) for y > 1. 

The longitudinal CF cr^^ itself, which also is shown in 
Fig. [TtI, has the same cusplike form with a gap described 
by (3.29), as in low di mensi ons. The linear q dependence 
in the denominator of ( 3.2£ ) says that, in spite of the gap, 
the correlation length of tr^^ near the surface is infinite 
at the anisotropic criticality. Actually, the correlation 
lengths near the surface are, in the ASM, the same as 
in the bulk. Indeed, above criticality cr^fj is a function 



of qz = ^Jk^ + cP] thus there are singularities in tr^^ at 
q — ±iKz, which cause the decrease of correlations of the 
type exp(— K^r) in the real space at large distances. 



VI. DISCUSSION 

In this paper, a comprehensive analysis of the behav- 
ior of the semi-infinite anisotropic spherical model at and 
above the ordinary phase transition is presented. The 
critical coupling of fluctuations, which usually necessi- 
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FIG. 14. Scaling function at zero wave vector, ^{x,0), in 
all dimensions. 



Extrapolation length 
in the semi-infinite ASM 
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FIG. 15. Extrapolation length vs lattice dimensionality. 



FIG. 16. Scaling function ^{x, y) of ( |1.18[ ) for d = 3. 
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FIG. 17. Scaling function F{y) of (4.22) for d = 2.5 



tales application of the e expansion or purely numerical 
methods, dies out in this model due to the infinite num- 
ber of spin components and makes in exactly solvable. 
On the other hand, the more important qualitative effects 
associated with Goldstone or quasi- Goldstone modes in 
weakly anisotropic magnetic systems are properly de- 
scribed by the ASM. The most important of these effects 
is the anisotropy-induced ordering in low dimensions. 
The ASM is superior with respect to the usual spherical 
model, which cannot incorporate anisotropy and yields 
unphysical results for spacially inhomogeneous systems 
because of the global spin constraint. On the other hand, 
the ASM is much better defined than its phenomeno- 
logical field-theoretical analog, the infinite-component 
model, and it can always be solved numerically. 

Unlike the renormalization-group (RG) approach, 
which is based on the expansion about the dimension 
d = 4 and becomes inefficient for low dimensions, the 
ASM describes the whole range 1 < < oo in a uniform 



way. The price for that is the rather complicated char- 
acter of the ASM system of equations in constrained ge- 
ometries, which makes application of numerical methods 
necessary. Nevertheless, there are a number of analytical 
solutions of the semi-infinite ASM in different limiting 
and particular cases. The most important of them are 
the isotropic-criticality solution of Bray and Moore for 
2 < d < 4, which was obtained above in an easier and 
more general way, the variations of the gap parameter 
Gin for c? < 2 and d > 4, and the slowly decaying gap 
tails of the correlation functions for g 3> k away from 
isotropic criticality. 

The gap parameter G„, or its deviation from the bulk 
value. Gin, plays a fundamental role in the theory of the 
ASM. The quantity — Gi„ is similar to the function V{z) 
used by Bray and Moore, and it also is proportional (and 
at criticality equal) to the inhomogeneity of the energy 
density, C/i„ [see ( 2.23| )]. The latter has been determined 
in Rcf. prt using RG and scaling arguments with the re- 
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suit Uln (x l/n*^^""^/" for 2 < c? < 4 at isotropic critical- 
ity, where v and a are the bulk correlation length and the 
heat capacity expo nents. In t he ASM 1/ = 2/{d — 2), as 
follows from (11^) an d I^AT\), an d a ^ {d ~ 4)/(d - 2), 
as follows from (|2.22 ) and ( 2.47 ). Thus the above for- 
mula reduces to Gi„ = C/i„ cx as was obtained 
by Bray and Moore. In these approaches, which con- 
sider n as a continuous variable, the inhomogeneous part 
of energy shows strong and unphysical divergence at the 
surface. Although it is clear that the continuous approx- 
imation is generally invalid near the surface, this strong 
singularity does not allow one to avoid the problem by 
replacing the surface region by some effective boundary 
condition, as can be done in the MFA. As a result, a 
numerical solution is principally ruled out for the semi- 
infinite field-theoretical 0(oo) model. In contrast, no 
such problems arise in the ASM, which is formulated on 
the lattice from the beginning. Moreover, continuous di- 
mensionalities (in the directions parallel to the surface) 
can be introduced into the ASM as well, while preserving 
the semi-infinite dimension discrete. The consideration 
in this paper shows that the nonscaling region near the 
surface, n 1, plays a very important role in the be- 
havior of the CFs in the asymptotic region far from the 
surface. So, the isotropic-criticality CFs are different for, 
say, d — 2.5 and d = 3.5, although they satisfy the same 
equation in the region n 3> 1. The difference between 
them stems completely from the region n ~ 1. 

The universality of the physical quantities in the ASM 
is different in different dimensionality ranges. For d > 4 
the gap parameter Gin is nonuniversal and decays as 
although the CFs have the universal mean- field 
form for n > 1. For 2 < < 4 both Gi„ and CFs 
are universal for n ^ 1 and nonuniversal for n ~ 1. 
For d < 2 the values of Gi„ are universal and decay as 
for n ^ 1. In contrast, the correlation functions 
are universal in the whole range of n. The reason for 
this strong universality and the ensuing strong scaling is 
that the (transverse ) CFs satisfy the constraint equation 
containing the integral over the wave vector q dominated 
by (7 ~ K <C 1 in low dimensions. 

The next steps in studying the inhomogeneous mag- 
netic systems with the help of the ASM should be (i) 
the solution of the semi-infinite problem below Tc and in 
field, (ii) inclusion of surface terms into the Hamiltonian, 
and (iii) numerical solution of the model in the film ge- 
ometry. The preliminary analytical investigation of the 
last problem can be found in Ref. [p4|. 
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